SEMI-CLASSICAL ORTHOGONAL POLYNOMIAL SYSTEMS ON NON-UNIFORM 
LATTICES, DEFORMATIONS OF THE ASKEY TABLE AND ANALOGS OF 

ISOMONODROMY 



N.S. WITTE 

Abstract. A B-semi-classical weight is one which satisfies a particular linear, first order homogeneous 
equation in a divided-difference operator D. It is known that the system of polynomials, orthogonal with 
respect to this weight, and the associated functions satisfy a linear, first order homogeneous matrix equation 
in the divided-difference operator termed the spectral equation. Attached to the spectral equation is a 
structure which constitutes a number of relations such as those arising from compatibility with the three- 
term recurrence relation. Here this structure is elucidated in the general case of quadratic lattices. The 
simplest examples of the B-semi-classical orthogonal polynomial systems are precisely those in the Askey 
table of hypergeometric and basic hypergeometric orthogonal polynomials. However within the D-semi- 
classical class it is entirely natural to define a generalisation of the Askey table weights which involve a 
deformation with respect to new deformation variables. We completely construct the analogous structures 
arising from such deformations and their relations with the other elements of the theory. As an example 
we treat the first non-trivial deformation of the Askey- Wilson orthogonal polynomial system defined by 
the q-quadratic divided-difference operator, the Askey- Wilson operator, and derive the coupled first order 
divided-difference equations characterising its evolution in the deformation variable. We show that this 
system is a member of a sequence of classical solutions to the q-Painleve system. 



Contents 

1. Background and Motivation 1 

2. Divided-Difference Calculus of the Quadratic Lattices 4 

3. Orthogonal Polynomial Systems on the Non-uniform Lattice 7 

3.1. General orthogonal polynomial systems 7 

3.2. g-Quadratic lattice 11 

4. Spectral Differences 13 

5. M = 2, L = case and the Askey- Wilson polynomials 22 

6. Deformation Differences 26 

7. M — 3, L — 1 case and deformation of the Askey- Wilson OPS 35 

7.1. Moments and Integrals 35 

7.2. Spectral Structure 38 

7.3. Deformation Structure 41 

7.4. Seed Solution 53 

8. Acknowledgments 55 
References 55 



1. Background and Motivation 

We propose a method for constructing systems of linear divided-difference equations which are analogs 
of isomonodromic linear differential equations and therefore "isomonodromic interpretations" or Lax pairs 
of the known g-Painleve and difference Painleve equations, such as those in the Sakai scheme [86]. What is 
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meant by an analog of a monodromy matrix turns out to be a connection matrix, appropriate to the class 
of lattices defining the divided-difference operators under consideration. In essence our method constructs a 
particular " isomonodromic analog" system from an orthogonal polynomial system, orthogonal with respect 
to a generalisation or deformation of a weight with discrete or countable support on a class of non-uniform 
quadratic lattices. We then deduce a number of linear divided-difference equations that this system satisfies 
and show that their pair- wise compatibility holds provided the coefficients of the linear system obey evolution 
equations of the difference or g-Painleve type. This is a very natural extension of the Fokas-Its-Kitaev 
construction [321 E31 E3] at the heart of Riemann-Hilbert techniques. 

Our method is independent of and distinct from other approaches which wc briefly recount here. The 
first studies to construct Lax pairs for the g-Painleve equations were those of Jimbo and Sakai |50j and 
Sakai [H7J HE] using the Birkhoff theory of linear g-diffcrcncc equations Y(qx) = A(x)Y(x) and imposing 
the condition that the connection matrix was independent of the zeros of det A(x). However this approach 
has not been extended beyond the D5 or cases (the latter case only found from a degeneration of 
the two- variable extension of the former). Another distinct approach which is founded upon the notion of 
the r- function of a rational d-connection is the Arinkin-Borodin theory [171 El 13] , which has been applied to 
the difference Painleve equations. Our approach is similar in spirit to that of Rains |85j . who has treated 
the master elliptic Painleve equation, in that an explicit construction is made of the solution to the linear 
problem which contains a multiple integral representation of the orthogonal rational function. However we 
will not make direct contact with this theory because it requires us to consider orthogonal rational functions 
on elliptic lattices, which is a generalisation beyond the class of lattices considered here. In addition we 
believe that the discrete Riemann-Hilbert approach, as formulated by Borodin in [16] and applied to two 
examples and extended to further cases in [18) . shares many features with the present study and one should 
be able to construct a Riemann-Hilbert formulation of our own theory. Recently Yamada has constructed Lax 
pairs for the q-Painleve equations for the systems with symmetries Eg , by a reduction from the 

elliptic form of the E^ Painleve equation [SSJ [§H] however no theoretical construction from first principles 
was proposed in the individual cases. 

The approach we propose here has been successfully employed for the isomonodromic systems and Painleve 
equations, see for example [66l EU [62l [13l [49l [19l EH ESI E31 [27l [94], and while most of the findings are 
reproductions of known ones it has lead to novel results hitherto not found using other treatments, such as 
the discrete Gamier systems |95j . An important feature of our approach is that it is strongly motivated by a 
probabilistic setting, namely that of the theory of random matrices and more generally determinantal point 
processes where the classical weights in the Askey table appear in the one-body factors of the eigenvalue 
probability density functions |25j . Some preliminary exploration of the program we propose here has already 
been initiated but not carried through to its logical conclusion, and we will delay citing this work until the 
body of our paper where it is directly relevant. However we should point out that we are most indebted to 
the pioneering work of Magnus [501 [S3] ■ 

The essential elements of our approach are the following. 

(A) The classification of special non- uniform lattices SNUL of quadratic type [60] , their associated divided- 
difference operators and ~M X and their rules of calculus, which applies to orthogonal or bi-orthogonal 
polynomial systems. Of significance is the fact that in general these lattices possess two fixed points which 
we denote Xl,xr- 

(B) The notion of a D-semi-classical weight w(x) which is characterised by an analog of the Pearson 
equation (see Definition I4.1[) [60l [63] 



where W{x), V(x) are polynomials in C[x]. 

(C) The orthogonal polynomial system (OPS) defined by such a weight on the SNUL Y n (x) £ C 2x2 
satisfies a three-term recurrence relation, which in our context is a particular Schlesinger transformation 
i!H)n+l (see [3T29]) 



(1.1) 



~D x w(x) 



2V{x) 
W(x) 



M x w(x), 



(1.2) 



Y n+ i{x) = K n (x)Y n (x). 



2 



(D) The semi-classical character implies, under fairly weak conditions, a spectral structure on the "isomon- 
odromic" system Y n (x) (60J [63] (see Proposition I4.4[) , namely that it satisfies the linear divided-difference 
equation 

(1.3) B x Y n (x) = A n {x)M x Y n (x). 

Here A n (x) is rational in x and the degrees of its numerator and denominator with respect to x are inde- 
pendent of n. 

(E) Parallel to the spectral structure is a deformation structure, whereby the weight and the system 
acquires a u dependence, constrained by another Pearson relation (see Definition 16. ip 

(1.4) V> u w(x; u) = 2 ^ X ' lU } m u w{x- u), 

R{x] u) 

where R(x] u), S(x; u) are polynomials in C[x]. 

(F) This also has a direct consequence for the "isomonodromic" system Y n (x; u), a deformation structure 
(see Proposition 16. 6[) and a second associated linear divided-difference equation 

(1.5) B u Y n (x; u) = B n {x; u)M u Y n (x; u), 

where B n (x;u) is also rational in x. 

(G) The compatibility relations implied by this over-determined "isomonodromic" system Y n (x;u) then 
lead to a number of conclusions - the deformation matrix B n (x; u) is expressible in terms of the spectral matrix 
A n (x;u) at neighbouring lattice nodes which we refer to as closure (see Proposition I7.8|) . and furthermore 
relations exist between components of the spectral matrix at two consecutive nodes on the u-lattice which, 
given a suitable parameterisation of this matrix by appropriate co-ordinates, is an recurrence relation on the 
deformation lattice. Our key results for such recurrence relations are given in Propositions 17.91 17.101 17.111 
and EE] 

(H) In our approach we derive a preservation property for the connection matrix defined as 

(1.6) P(x;u) :={Y R {x;u))- 1 Y L {x;u), 

where Yl^x; u) are local fundamental solutions of the spectral equation about xl,r- Specifically this implies 

(1.7) B x P(x;u)=0, H u P(x;u)=0, 

as a consequence of our assumptions. 

In our approach we can sidestep a number of issues, which are important to be sure, but don't affect 
the outcome. So we postpone deeper considerations of an analytic or algebraic theory of linear systems of 
divided-difference operators on these lattices for subsequent studies, however our present work will provide 
concrete illustrative examples for such an investigation. We are referring to, for example, issues of a Galois 
theory for D-difference equations, a Birkhoff theory for the local character of the solutions to systems of 
linear first order D-diffcrcnce equations and analogs of monodromy for D-difference equations. 

Our approach poses the question concerning a correspondence between the system of hypergeometric 
and basic hypergeometric orthogonal polynomial systems generalising the classical systems, known as the 
Askey table, as represented in reference work [56] . and the most complete system of elliptic, q-difference 
and difference analogs of the Painleve equations, the Sakai scheme [86]. A correspondence between the 
Sakai scheme and the Askey table would explain the occurrence of many features that have been discovered 
recently such as the appearance of basic hypergeometric functions in their classical solutions. We have 
already referred to the correspondence between the classical orthogonal polynomial systems (OPS) Hermite, 
Laguerre and Jacobi and the classical solutions to Painleve IV, V and VI respectively. Amongst the difference 
and g-difference OPS the current evidence for such a correspondence can be summarised in the Table [T] 

The layout of this work loosely follows the plan given above. 

In Section [2] we describe the classification of the quadratic lattices and their divided-difference calculus, 
with special emphasis on the example of the master class, the g-quadratic lattice. Section[3]is devoted to the 
reformulation of orthogonal polynomial system theory necessary for systems with a weight having support on 
a general quadratic lattice, and some analytical results for the g-quadratic lattice. The spectral structures 
are laid out in Section^] for a general quadratic lattice, starting with the definition of a D-semi-classical 
weight and developing the consequences of this for the orthogonal polynomial system. We also describe 
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Base OPS 


Integrable system 


Reference 


little g-Jacobi 


Dg', full g-Pvi 


m 


Pastro [80] 


full g-Pvi 


M 


little g-Jacobi 


L>5 , special q — Pyi 


m 


g-Krawtchouk 


L>5 , special q — Pyi 


us 


q-Charlier 


degenerate q — P V i 


[13 


g-Freud 


degenerate q — Py 


m 


Meixner, Krawtchouk 


d-Py 


m 


Charlier 


d-Pw 


m 



Table 1. The Askey table-Sakai scheme correspondence ranked according to the degenera- 
tion pathway from the master case. The afhne Weyl group refers to the symmetry group of 
Backlund transformations for the integrable system. Only the discrete lattices are included 
here. 



the compatibility relations of this structure with the three-term recurrence relations, and this leads to a 
generalisation of the Laguerre-Freud equations. An explicit example of the forgoing theory is the lowest 
case, and from our point of view the trivial case, of the q-quadratic lattice which is denoted by the label 
M = 2. This is dealt with in Section [5] We recover the Askey- Wilson OPS and demonstrate that every 
aspect of this system can be derived in an efficient manner using the theoretical tools developed in the 
previous section. 

In Section [6] we introduce an (or many) auxiliary variable and demand that the weight satisfies an analo- 
gous D-semi-classical relation to the spectral relation, on a general lattice not necessarily the same as that for 
the spectral variable. This constitutes a deformation of the weight which leaves it having the same structure 
as the original. From this we develop a parallel theoretical analysis of the orthogonal polynomial system 
with respect to this deformation. We now have two compatibility relations with this deformation structure 
- one with the three-term recurrence and other with the spectral structure. For the present purposes this 
completes our task for the general theory. 

In Section [7] we treat a natural deformation of the Askey- Wilson weight, but not by any means the 
only one, to the q-quadratic lattice. Here we construct parameterisations of the spectral and deformation 
matrices, the closure relations between the two structures, determine the co-ordinate transformations to a 
set of variables and derive the evolution equations for this system (see Propositions 17.91 l7.10|7.lH and l7T2|l . 
Finally we offer evidence that our system is a classical solution of the £y q-Painleve system. 

2. Divided-Difference Calculus of the Quadratic Lattices 

Divided-difference operators associated with the special non-uniform lattices have appeared in many 
studies of orthogonal polynomials of a discrete variable. For example see the early studies by Hahn [SHI 
l35l l37l [38] . the foundational work by Askey and Wilson [4] and the monograph of Nikiforov, Suslov and 
Uvarov [78] . The employment of divided-difference operators such as the Askey- Wilson and Wilson operators 
has been common place in studies of the Askey table of hypergeometric orthogonal polynomial systems too 
numerous to cite here (see the monographs by Ismail [43] and Lesky [59]). 

Let II„ [x] denote the linear space of polynomials in x over C with degree at most n € Z>o . In pioneering 
investigations Magnus |60[ 163) provided a geometrical understanding of these lattices and their divided- 
difference operators which we relate here briefly. If we define the divided- difference operator (DDO) 

by 



(2.1) 



then a simple consequence of the condition that : n„[x] — > n n _i[x] for all n 6 N is that l±(x) are the 
two j/-roots of the quadratic equation 

(2.2) Ay 2 + 2Bxy + Cx 2 + TDy + 2£ x + T = 0. 
The functions l±(x) satisfy 

r~ ~s / \ / \ „<Bx + 2? 

(2.3) t + (x) + t_(x) =-2— — , 

, , , , Cx 2 + 2£x + T 

(2.4) H (x>_(x) = , 

and their inverse functions i^ 1 are defined by t± (t-j-(x)) = x. For a given y- value the quadratic (|2.2[) 
defines two x-roots if C ^ 0, which are consecutive points on the x-lattice, x s := x(s),x s +i :— x(s + 1) 
parameterised by the variable s and therefore defines a map x s i— > x s+ i. Conversely for a given x- value 
the quadratic defines two y-roots if .4 =/= which are consecutive points on a dwa/ lattice, the y-lattice, 
y s := y(s) = L^(x(s)),y s+ x :— y(s + 1) = l + (x(s)) which generally is distinct from the x-lattice. We 
will employ an operator notation for the mappings from points on the direct lattice to the dual lattice 
E±f(x) :=/Mx)) so that 

(2.5) B x f(x) := /M*)) -/(*-(*)) = Etf-E-f 

i+{x)-L-{x) E+X-ExX 

for arbitrary functions /(x). The inverse functions l^}(x) define operators (E^)^ 1 which map points on the 
dual lattice to the direct lattice and also an adjoint to the divided-difference operator 

t + (x) - t_ (x) (EJ) X X - (i/a, ) X 

The composite operators E x := (E~)~ 1 E+ and -E^T 1 = (E+)~ 1 E~ map between consecutive points on the 
direct lattice. 

However in the situation of a symmetric quadratic A = C and T> = £ , which entails no loss of generality, 
then we have (i?+) _1 = E~ and (K,?) -1 = Consequently there is no distinction between the divided- 
difference operator and its adjoint, and hereafter we adopt this simplification. For example one useful 
consequence of this choice is that the ratio 

(x - E+u)(x -E~u)_ Ax 2 + 2Bxu + Cu 2 + 2Vx + 2£u + T _ 
^ 2 ' 7 ' (E-x-u)(E+x-u) ~ Au 2 + 2Bux + Cx 2 + 2Vu + 2£x + T ~ 

for all x, u. A companion operator to the divided-difference operator D x is the mean or average operator M x 
defined by 

(2.8) M x /(x) = i[/Mx)) + /( t _(x))], 

so that the property M x : II„[x] — > II n [x] is ensured by the condition we imposed upon D x . 

Definition 2.1. Henceforth we use the shorthand for the difference in consecutive y points Ay(x) := 
t+(x) — t_(x). We introduce the notion of "fixed points", whereby Aj/ 2 (xj?) = 0, or L + (xp) — L-(xf) and 
are given as the roots of (B 2 - AC)x 2 F + 2(BV - A£)x F + V 2 - AT = 0. 

Assuming AC ^ one can classify these non-uniform quadratic lattices (or SNUL, special non-uniform 
lattices) according to two parameters - the discriminant B 2 — AC and 

(A B 

(2.9) 6 = det \B C £ 

\V £ T/ 

or AQ — [B 2 — AC)(T> 2 — AT) — (BT) — A£) 2 . There are four primary classes - the linear lattice, the linear 
g-lattice, the quadratic lattice and the g-quadratic lattice, which are given in Table [21 The g-quadratic 
lattice, in its general non-symmetrical form, is the most general case and the other lattices can be found 
from this by limiting processes. 

This classification of lattices for polynomial systems can be extended to rational function systems [89l l90l 
I64j and in this case one has a lattice characterised by a bi-quadratic relation and parameterised by elliptic 
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B 2 -AC 


e 


conic 


lattice 


canonical DDO 


Notes 








parallel lines 


linear 


forward difference 




> 





intersecting lines 


g-linear 


g-difference 







< 


parabola 


quadratic 


Wilson 




> 


< 


hyperbola 


g-quadratic 


Askey- Wilson 


q real 


< 


< 


ellipse 


g-quadratic 


Askey- Wilson 


M = i 



Table 2. The non- uniform lattices of quadratic type. 



functions [33] (see also Section 15.10 of [H]). However for the purposes of studying the Askey table we will 
not pursue this direction. 

For the quadratic class of lattices the parameterisation on s can be made explicit through the trigonomet- 
ric/hyperbolic functions or their degenerations so we can employ a parameterisation such that t_(a;(s)) = 
y(s) = x(s — 1/2) and l+(x(s)) = y(s + 1) = x(s + 1/2). We denote the totality of lattice points 
by G[x] := {x(s) : s € Z} with the point x(0) — x as the basal point, and of the dual lattice by 
G[x] := {x(s) : s G 

Having established the basic properties of the divided-difference operators we can deduce key elements of 
their calculus. A consequence of the general definition of the divided-difference operators are the following 
identities 

(i) the product or Leibniz formulae 

(2.10) B x fg = B x fM x g + M x fB x g, 

(2.11) M x fg = M x fM x g + \Ay 2 B x fB x g, 

(ii) the inverse formulae 

1= . D .i %/ 



(2.12) 



E x fE x f 



f E+fE x f' 

(iii) and the commutativity formulae 

(2.13) M U M X = M X M U , M U B X = B X M U , B U B X = B X B U . 

The other side of our divided-difference calculus concerns the definition and properties of analogs to 
integrals. We define the O-Integral of a function defined on the x-lattice / : G[x] — > C with basal point xq 
by the Riemann sum over the lattice points 



(2.14) 



f Bxf(x) := 5>+(*) - t-(s))f(x(s)) = ^s)f(x s ), 



where the sum is either a finite subset of Z, namely {0, . . . ,91}, or Z>o, This definition reduces to the 
usual definition of the difference integral and the Thomae- Jackson g-integrals in the canonical forms of the 
linear and q-linear lattices respectively. A number of properties flow from this definition - 

(i) an analog of the fundamental theorem of calculus 



(2.15) 



DxB x f(x) = f(E+x m ) - f(E-x ), 



XQ<X s <X-.-yx 



(ii) an analog of integration by parts for two functions f(x),g(x) 
(2.16) 

f Bxf(x)B x g(x) = - f B(E+x)B x f(E+x)g(E+x) + f(E+ 2 x < n)g(Etx < n)-f(x )g(E x xo), 



(iii) and the parameterisation "invariance" property 



(2.17) / Bxf(x) = B(E+ 2 x)f(E+ 2 x) + Ay(x Q )f(x )~Ay(E+ 2 x <n )f(E+ 2 x m ), 

J XqKx s KX(JI J Xq<X s <X?)i 



(2.18) = / B(E-'x)f(E-'x)-Ay(E-'xo)f(E-'x ) + Ay(x<n)f(xn). 

J Xq <X s <Xyx 

We will apply our theory to the case of the q-quadratic lattice and the Askey- Wilson divided-difference 
calculus, and in order to simplify the description and to conform to convention we will employ the canonical, 
that is to say the centred and symmetrised forms of the lattice and the divided-difference operators. Let us 
define the base q = exp(2ir]) although we will not restrict ourselves to q-domains such as < $i(q) < 1 except 
to avoid special degenerate cases and to ensure convergence. Consider the projection map from the unit 
circle z = e l8 , 9 e [— tt, it) onto [—1, 1] by x = \{z + z^ 1 ) = cos 9 e [—1, 1]. We denote the unit circle by T 
and the unit open disc by D. The inverse of the projection map defines a two-sheeted Riemann surface, one 
of which corresponds to the interior of the unit circle, and the other to the exterior. Thus we take a;-planc 
to be cut along [—1, 1] and will usually give results for the second sheet i.e. when \x\ — > oo as z — > oo. In 
the symmetrised and canonical form of the lattice we have A = C, arbitrary and non-zero, B = — cos r\A. 
V = £ = 0, T = -sin 2 77.4, and 9 = 2s?y. Define the shift operators E± by E±f(x) = /(^^z + q-^-^z" 1 }) 
and set y± = E^x. This implies that 

(2.19) V++V- = (q 1/2 + q~ 1/2 )x = 2cosr) x, 

(2.20) Ay := y+ - y_ = \{q 1/2 - <T 1/2 )(z - z" 1 ) = -2 sin 77 sin 0, 

(2.21) Ay 2 = (q 1 / 2 ~q- 1 / 2 ) 2 (x 2 -l), 

(2.22) y + y_ = x 2 + {(q 1 ' 2 - q' 1 ' 2 ) 2 = x 2 - sin 2 r,. 
The Askey- Wilson divided-difference operators are defined as 

, , , , fiUq^z + q-^z- 1 ])- fa[q- l / 2 z + q l ' 2 z- 1 ]) 

(2 - 23) ^ = ^ WI>-%')(z-*-i) " 

M x f(x) = \ [f^z + q-^z-l]) + f{\[q-^ 2 z + q^z' 1 }) 

There is an explicit parameterisation of the q-quadratic lattice 

(2.24) x{s) = \{q s + q- s ) = cos(2r/s), 

(2.25) y±(s) = \{q s±1 ' 2 + q-^ 1 ' 2 ) = cos(2r,[ S ± |]). 

Here the direct lattice is G[x — j{a + a -1 )] = {^(q r / 2 a + q~ r / 2 a~ 1 ) : r e 2Z}. For \q\ — 1 and r/ not a 
rational multiple of 7r then the lattice densely fills the interval [—1, 1]. If 7/ is a rational multiple of 7r then 
one has the root of unity case q N = 1 and a finite lattice. In the generic case we will assume that we are 
dealing with functions f(x) in the class where 

(2.26) J Hxf{x) = S -^- J dxf(x), 

is applicable. The reader should note that we will not distinguish a function of x, f{x), from the function of 
z, f{z) = f(^[z + z^ 1 ]) as done by some authors, and it should be clear from the context which is meant. 

3. Orthogonal Polynomial Systems on the Non-uniform Lattice 

3.1. General orthogonal polynomial systems. Our study requires the revision of a number of standard 
results in orthogonal polynomial theory EH1 S3] so we recount our formulation. Let {l n (x; a)}^ =Q be a 
polynomial basis of L 2 (w(x)U>x, G), where is of exact degree n and the support is G — {E+ k x : k g 2Z} 
or if finite G — {xo, ■ ■ ■ , x<ji} and a denotes the set of parameters characterising the lattice. The appropriate 
canonical basis is dependent on the lattice type through the general requirements that 

(i) l n is of precise degree n so that l n (x; a) = g n (a)x n + 0(x n ~ 1 ) with g„ ^ 0, 
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(ii) is an exact lowering operator in this basis 

(3.1) ®> x ln(x; a) = Cn{a)l n -x{x] a'), 

where c„ is constant with respect to x and the transformed parameter set a 1 is related to the original 
a depending on the lattice type. 

We will also require the linearisation formula xl n (x; a) — d n (a)l n +i(x; a) + e n (a)l n (x; a). A general solution 
to the two requirements above is the following product expression 

n-l 

(3.2) l n (x; a) = g n (a) [J [x - (E+) 2k x(a)}, 

k=0 

where the basal point x{a) is parameterised by a. We note that for some lattices the limit lin^^oo l n (x; a) 
exists and in this case we denote it by ioo(^i a )- For the classes of quadratic non- uniform lattices the basis 
choices are tabulated in Table OH 



DDO 


Lattice Type 


Basis l r 


Notes 


d 

dx 


continuous 


x r 


c r = r 


A x 


linear 


^ -%-")- ^". + T \ 

fc=0 v ; 


c r = r 


D q 


q- linear 


(ax;q) r = l\(l-aqx)= iaqrx;qU 


1 - aq r 

c r — 

q - 1 

a' = qa 


W 


quadratic 


TT \t I (h I n\ 2 ] - + a ~ *^) r ( r + a + ^) 

lU-+( fc + a ) \ T{a-i^)Y{a + i^) 


c r = T 
a! = a + i 




g-quadratic 


{az,az ;q) T = 

(aq r z,aq r z ,g)oo 


<? - 1 
c r ~ —2a 

9-1 
a' = q^a 

g r = {-2a) r q^ r[r - 1] 

d r = 

2aq r 

e r = \(aq r + a~ 1 q~ r ) 



Table 3. Canonical Bases for the non- uniform lattices of quadratic type. 



Consider the general orthogonal polynomial system {p n (x)}^ =0 defined by the orthogonality relations 

(3.3) / Dx w(x)p n (x) l m (x; b) = \ ° ® ~ 171 <n ; n>0, 

Jg [hnK°) m = n 

with G denoting the support of the weight w{x). Our system of orthogonal polynomials and their associated 
functions (to be defined in (|3.21[1 ) have a distinguished singular point at x = oo and possess expansions about 
this point which can characterise solutions uniquely. This is related to the fact that orthogonal polynomials 
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are the denominators of single point Pade approximants and that point is conventionally set at x — oo. We 
give special notation for the coefficients of x n and a;™ -1 in p n (x), 

Pn(x) = l n X n + -Jn^X^ 1 + . . 



n > 0. 



(3.4) 

The corresponding monic polynomials are then ir n (x) = 7~ 1 p„(a;) given that n > 0. A consequence of the 
orthogonality relation is the three term recurrence relation 



(3.5) 



a n+1 p n+ i(x) = (x- b n )p n (x) - a n p n -i(x), n > 0, 



and we consider the set of orthogonal polynomials with initial values p_i = and po = 7o- The three term 
recurrence coefficients are related to the leading and sub-leading polynomial coefficients by [HH [3D] 

_ 7«-i , _ 7n,i _ Jn+hl 

7n In 7n+l 



(3.6) 



n > 1. 



The initial values of the recurrence coefficients are 
(3.7) 



oo = , 7oa = 0. 



7i 

and where ao is not fixed by the initial polynomials but rather by the initial associated functions (see after 

(E22D. 

The orthogonality relation (|3.3p is derived from the linear functional on the space of polynomials £ : p € 
II i—» C and we employ our basis polynomials as an expansion basis although not necessarily with the same 
parameter as above 



(3.8) 



Pn(x) = >^ Cn,fc(a)/fc(a;; a) n>0. 



k=0 



Consequently we define the moments {mj t k}j,k=o,i,...,oa of the weight as the action of this functional on 
products of the basis polynomials, defined as 



(3.9) 



mj t k(b,a) 



Dx w(x) lj(x;b)lk(x; a), j,k>0. 



Central objects in our theory are the Moment determinants 

(3.10) A„ := det[mj-,fc]j,fc=o,...,n-i) n>l, A := 

and 

mo,o • • • m o,j-i [] moj+i • • • mo t n 



1, 



(3.11) S n>j := det 



n > 1, j = 0, . . . n — 1. 



"In— l.Ti 



defined in terms of the moments above. Obviously A n = S„. n and we set Eo,o := 0. The expansion 
coefficients are given in terms of these determinants 



(3.12) 

ft follows from 



that 



E A. 
c n . 3 (a) = (-)™ + ^„(o)- r ^ iL , c„.„(a) = K(b)—^— 



IX w(x)[p n (x)} 2 = Cn n (b)hn(b), n > 0, 



and thus for p n {x) to be normalised as well as orthogonal we set c njn (a)h n (a) = 1. We have moment 
determinant representations of the polynomials 

/ m ,o • • • mo,j ■ ■ ■ m ,„ \ 



(3.13) 



Pn{x) = Cn '"^ det 



m n -i t o 
V ^o 



m n -i,j 
U 



m n -l,n 
in. 



n>0. 



The three-term recurrence coefficients are related to these determinants 



(3.14) 
(3.15) 
(3.16) 



a 2 n = dn-i(a)dn-i(b) 



K = e n (a) + d n (a) " +1 ' n - d n -i(a) 



n>l, 

z^n.n— 1 



n > 0, 



iri+l 



li = 9n{a)g„(b) 



A 



n > 0, 



Tt + l 



where each coefficient is independent of the choices of a, b as can be easily verified from their determinantal 
definitions given previously. 

Another set of polynomial solutions to the three term recurrence relation are the associated polynomials 
{p£\x)}%L , defined by 



Pn-l( X ) 



I sPn(y)-Pn(x) 

y w(y) , ri > 0. 



y-x 



(3.17) 

In particular these polynomials satisfy 

(3-18) a n+ ip£\x) = (x- b n )p%-i(?) - a n p { ^L 2 {x), 

with the initial conditions p^\(x) — 0, p$ (x) = ?no,o7i. Note the shift by one decrement in comparison to 
the three-term recurrence (|3.5[) for the polynomials {p n {x)}^ - We also need the definition of the Stieltjes 
function 



x-y 



(3.19) f{x) = 

which is a moment generating function in the following sense 

;a) v-^ m ,„(a) 



x i G, 



(3.20) 



f(x) 



E 



x 4. G, x — > oo, 



(x;a) 4^ d n (a)l n +i(x; a) ' 

and splits into two parts - one part being a series with inverse basis polynomials and a remainder, which 
may be absent for some lattices. We define non-polynomial associated functions or functions of the second 
kind {q n (x)}%L by 

(3.21) q n {x) :=f{x)p n {x)-p^ ) _ l {x), n > 0, 
which also satisfy the three term recurrence relation (|3.5[) . namely 

(3.22) a n+ iq n+ \(x) = (x - b n )q n (x) - a n q n -i(x), n > 0, 

subject to the initial values q-\(x) — l/ao7o, qo(x) = 70/(2;)- The initial value of ao is irrelevant and 
therefore arbitrary in so far as the polynomials are concerned, however many relations for the whole system 
extend from n > 1 to include n — if we allow this to be finite, non-zero and satisfying the above initial 
condition. The associated functions also have a determinantal representation 

/ m , • • • m , 3 ■ ■ ■ ™o,« \ 

.(a) 



(3.23) 

where 
(3.24) 



q n {x) 



det 



V fo 



m n -x t j 

Si 



'it'n — l,rt 
fn J 



n > 0, 



* - V 



fj(x;a) := 

Likewise this function has an expansion analogous to p.20j) 

foo,j(x;b) ^ m J> (&, a 



, j > 1, /o = .9o/- 



(3.25) 



fj(x;b) 



^00 ft) 



E 



d„(a)/ n+ i(a;;a)' 
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x £ G, x 



The polynomials and their associated functions satisfy the Casoratian relation 



(3.26) 



1 



p n {x)q n -i{x) - p n -i(x)q n {x) = — , n > 



Central to our analysis is a composite of polynomial and non-polynomial solutions of p.5[) . the 2x2 
matrix variable 



(3.27) 



Y n (x) = 



Pn(x) 



q n (x) 
w(x) 



n > 0. 



We will refer to this as the orthogonal polynomial system (OPS). From (j3.26p we note that 

1 



(3.28) 



det Y n (x) 



n > 0. 



a n w(x) 

The three term recurrence is then recast as the matrix equation 
(3.29) Y n+1 (x) =K n Y n {x), n > 0, 

with the recurrence matrix given by 



(3.30) 



K n (x) 



1 I x bn dji 




det K„ 



n > 0. 



Ctn+l 



a n +i \ a n+l 

A well known consequence of (|3.5[) are the Christoffel-Darboux summation formulae 

n-l 



(3.31) 
(3.32) 
(3.33) 



j=0 
n-l 

J2q j {x)p j {y) 

n-l 
3=0 



[Pn{x)p n -l{y) ~ Pn-l(x)Pn{y)\ . n 

a n , n > 0, 



x-y 

[q n {x)p n - 1 (y) - q n -i(x)p n {y)] 



1 



n > 0, 



a; — y x — y 

[q n {x)q n -i(y) - q n -i{x)q n (y)} f(x) - f{y) 



n > 0. 



x-y 



x-y 



If one is only interested in the leading orders of the large x expansion rather than a systematic expansion 
then it is convenient to employ an expansion in monomials rather than in basis functions. Extending (|3.4[) 
we have expansions about the fixed singularity at x = oo 



(3.34) 



Pn{x) = 7n 



/n-l 



n-l 



■,i=0 / \Q<i<j<n i=l 



valid for n > 1, while for the associated functions 



(3.35) 



9n(z) = 7n' 



n+l 



-n— 2 



\i=0 



E & a+E' 



-n— 3 



+ 0(a;-"- 4 ) 



valid for n > 0. 



3.2. g-Quadratic lattice. In Sections [5] and [7] we intend to apply our theory to the g-quadratic lattice and 
will draw upon numerous properties of the corresponding basis, which we discuss here. Firstly we recall the 
analytic continuation of (j> n (x; a) which is 



(3.36) 



4> r (x;a) 



(az,az 1 ;q) c 



(aq r z,aq r z 1 \q) Q o' 

for all r G C but subject to |g| < 1. We will employ the shorthand notation (az ±:L ; q)^ = (az, az~ x \ q)oo. 
Implicit in the above formula is the elliptic-like function (/^(a^a) which has meaning for all a,x € C for 

li 



\q\ < 1. The action of the divided-difference operators are 

(3.37) ® x (/> r (x; a) = 2a ^— y <t>r-i(x; q 1/2 a), 

(3.38) M*^.(a;;a) = |(1 + q~ r )Mx; q 1/2 a) + |(1 - g~ r )(l - oV r "%-i(iti Q 1/2 a), 

(3.39) IMoofo a) = -^-<Mz; g 1/2 a), 

9-1 

(3.40) Mx^fo a) = ^(x; <T 1/2 a). 

The expansion theorem of Ismail [40] [41] states that 
Proposition 3.1. Let p be a polynomial of degree n, then 

n 

(3.41) p{x) = ^2 Pk<t>k (x; a) 
for any a € C where 

(3-42) p k = {q - 1)k g- k ( k -W(B*p)(x k ) 

{2a) k {q; q) k 

with Xk = \(q k / 2 a + q^ k ^ 2 a^ 1 ). 

In these works, see (2.2) and the proof of Theorem 1.1 in [3D], we also have the change of base formula 

" n ( b \ k 

(3.43) 4>n{x\ b) = ^2 [™] {abq k , b/a;q) n ^ k y-J tj> k (x;a), 

where we use the standard definition of the g-binomial coefficient. This result allows us to derive the following 
linearisation formula. 

Lemma 3.1. The product of two basis polynomials with the same base a has the following expansion in 
terms of the same basis 

(3.44) Mx;a)Mx;a) = (-) k+l q- kl £ (-) m (a 2 q m ; q ) k+l _ m — ^,q)i(q\q)k m (x;a), 

m= ^ (k l) q)k+i-m(q; q) m -k{q; q) m -l 

for all k, I e N. 

Lastly we have the Cauchy expansion formula of Ismail and Stanton [47j . 
Proposition 3.2. The Cauchy kernel has the expansion 
(3 45) 1 = 1 0°o(x;a) 2a y^ (j>n(x;a) ^ n 

for all y such that 1/ / 1 and </>oo(j/;a) 7^ 0. The expansion also holds for y — j/o with </>oo(yo;a) = but 
yo 7^ x in the sense that the left-hand side (j/o — a;) -1 equals the limit of the right-hand side as y — > j/o- 

This implies the following expansion of the Stieltjes function. 

Corollary 3.1. The Stieltjes function has the following expansion as x — > 00 with x 7^ x k (a), with k € 2Z 
(i.e. 4>oo(x; a) 7^ 0) 

(3-46) f(x) = - 2a jr - g " . m , n (a), 

9oo{x;a) ^ (p n+1 {x]a) 

where 

(3.47) f 00 (x)= ( Byw(jj)^ t!r " ] 



g x ~v 
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Conforming with standard notation |31) we define the basic hypergeometric function r +ifr by the series 



(3.48) r+Wr 



a±, a 2 , a r +\ 

bu b 2 , W ■ A ' Z 



E 



(qi, a 2 , . . . , a r +i; q) n ^ n 
(q,bi,b 2 ,...,b r ;q)n 



which is convergent for \z\ < 1. The very-well-poised basic hypergeometric function r+ \W r is a specialisation 
of the above 



ai, q^/ai, -q^/ai, 04, a r+ i 



• Q % 

at, -y/ai, qai/cLA, qai/a r+ i ' 



(3.49) r+ i W r (ai, a 4 , a 5 , . . . , a r+i ; z) = r+1 tp r 



so that qai — a 2 b\ = a^b 2 = ■ ■ ■ = a r+ \b r and a 2 — q^/ai, 03 = — q^/ai. The r+i^fir or r +iW r functions may 
also be balanced, whereby z — q and Ilj=i bj — q Y\j=i a j- 

4. Spectral Differences 

In this section we lay out the structures of the spectral divided-difference operator in the context of the 
orthogonal polynomial system for a general lattice type. Our starting point is the notion of the D-semi- 
classical weight, as given by the following definition of Magnus [60 . 

Definition 4.1 ( 6Q\). Let the D- semi- classical weight satisfy 

(4.1) WB x w = 2VM x w, 

or equivalently 

(4 2) My^^ W + AyV 

for W(x), V(x) irreducible polynomials, which we will call spectral data polynomials. Furthermore we assume 
W ± AyV for all x E G. For minimal degrees of W, V this is the analog of the Pearson equation. 

Remark 4.1. On the finite lattice x G {xq, . . . ,x<yi} we naturally require w(x) ^ 0, however we will impose 
upper and lower terminating conditions 

(4.3) w(E+ 2 x m ) = w(E- 2 x ) = 0, 

respectively. These are the conditions analogous to Eqs. (2.3.2) and (3.3.3) in Nikiforov, Suslov and Uvarov 
[78] . To be consistent with (|4.2[) this implies 



(4.4) (W + AyV)(E+x<n) = (W - AyV)(E- Xo ) = 0. 

Remark 4.2. In a series of works Suslov and collaborators [5 , 91, 6, 84, 82) have studied the Pearson equation 
for all of the lattices admissible in the classification, and sought solutions for the weight functions given 
suitable polynomials for W, V. However they limited their choices to examples of minimal degree for W, V 
which made contact with those weights found in the Askey table. In terms of our own variables those found 
in [9T] are given by 

(4.5) a = E-(W-AyV), T = K(W + AyV)( X ) -E~(W - AyV)( X ) 

From (|3.20[) we recognise f(x) as a moment generating function and a key element in our theory are the 
systems of linear divided-difference equations satisfied by the moments. Before we state such systems we 
need to note the following result. 

Lemma 4.1. Let l^{x;a) be a canonical basis polynomial. For any k £ Z> and a £ C the weight satisfies 
the integral equation 

1 w i x) 

(4.6) / Bx-y,{(l k [W + AyV})(E+x)-(l k [W-AyV})(E-x)}=0. 

Jx <x<x m ^y\ x ) 
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Proof. We will prove this statement accounting for the boundary terms explicitly. We establish two prelim- 
inary identities first, namely 



(4.7) 
(4.8) 



3>x w(x) 



}xw(x) 



1 



& _ <(ww+asv1) __^_ 



Xq <X<Xyx 



\Ay(E+AyY 



[W + AyV]) 



Ex Ay 



E-(l k [W-AyV])j=0. 
E- (l k [W - AyV})\ = 0, 



with l k :— l k (x; a). Taking the first of these we compute 



= 



T>xl k {E+x)(WB x w - 2VM x w)(E+x) 



Xq <X<Xy 



Ay(x s )l k (E+x s )W(Ei 



w(E+ 2 x s ) - w(x s ) 



Q<s<<Xt 



E 



Ay(E+x s ) 

Ay{x s )l k {E+x s )V{E+x s ) [w{E+ 2 x s ) + w(x s )] , 

o<s<<n 

^y(E- 2 x s ) 



E 



Ay{x s 



Ay(E^xo) 



l k (E+x s )W(E+x s )w(x s ) 



l k (E-x )W(E-x )w(x ) 



Ay(Ex x<y{) Ay{ExX Q ) 
Ay{E- 2 x s )l k {E-x s )V{E-x s )w{x s )- ^ Ay{x s )l k {E+x s )V{E+x s )w{x s ) 



0<s<9t 



Ay(x m )l k (E+xvi)V(E+xvi)w(E+ z xsjx) + Ay{E- 2 x a )l k (E^ xo)V(E~ x )w(x ), 



^ Ay(x s )w(a 



0<s<9t 



- I 



w 
Ay 



V 



(E+x s 



&y(E~ 2 x s ) 



W 



Ay 



- V 



(E-x.) 



A ?l X + n) J k (E+x n )w{E+*x n )[W - AyV](E+x m ) 



Ay(E+x m ) 



&y(E- 2 x ) 

Ay(E x x ) 



l k (E-x )w(x )[W - AyV](E-x ) 



where the boundary terms arise from using the integration by parts and the change of variables x s i-> E^ 2 x s 
formulae. Both boundary terms are in fact zero due to the upper and lower terminating conditions. The 
second formula (|4.8p is just a variation of the first with the integrand evaluated at E~x s instead. Now for 
a general quadratic lattice 

Ay + E^Ay _ 

E±Ay A' 

so we simply add (|4.7j) and (|4.8[) and note that a common, constant factor appears in the integrand of 
precisely this form. Then (|4.6|) immediately follows. □ 



Remark 4.3. An equivalent system for general lattices has been derived by Suslov [91] in the special case of 
Jfc = of I 



We also need to express W ± AyV in terms of canonical basis polynomials and formulate the following 
definitions for the coefficients Kjv,z(a), <Jjv,i(°) 



N 



(4.9) 



(4.10) 



E+(W + AyV)(x) + E-(W - AyV)(x) =: 2 £ n NM {a)l k {x- a), 

N-l 

E+(W + AyV)(x)-E-(W-AyV)(x) =:AyJ2S Nik (a)l k (x;a). 



k=0 



Here N is the cutoff determined by the degrees of W, V and the lattice type. 
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Corollary 4.1. Let us assume that b' ^ a, where b,b' are related as a, a' are by \3.1\) . For all k £ Z>o the 

moments m kj i(b',a) are characterised by the linear, homogeneous recurrence relations 



(4.11) 



d k (b)c k+1 (b) + ^d k ^(b')c k (b) 



N-l 



^ S N: i(a)m ki i(b',a) 



1=0 



N 



1=0 



B V 

kn,i{cl) + [ e k (b) + — e k -i(b') + ) S N ,i{a) 



m k -ij(b',a) = 0. 



For k — the last term is absent and therefore b, b' do not appear. 



Proof. We start with the integral formula (|4.6p and employ the resolutions of E^(W + AyV) and E~(W — 
AyV) as given by (|4.9I) and (|4.10p . In addition we require expressions for H) x l k (x; b) and M x l k (x; b), and use 
(I3.ip for the former and 
(4.12) 



B 



A 



J k (x;b) = d k (b)c k+1 (b) + -d k -x{b')c k {b) l k {x;b') + c k (b) e k {b) + -re fc _i(6') + -r Jfc-iM'), 



B 



for the latter. 



□ 



If a specialisation of the parameters b',a is made to effect some cancellation with corresponding factors 
in the weight then the moments are given by an integral whose integrand has the same structure and (14. lip 
is now a linear divided-difference equation with respect to the internal parameters of the weight. 

Now we continue to develop the consequences of the ID-semi-classical weight on the orthogonal polynomial 
system, and in particular for the Stieltjes function. The following result will be crucial for our approach. 

Proposition 4.1. Given a D- semi- classical class of an orthogonal polynomial system fulfilling \4- 1\ ) and the 

conditions therein then the Stieltjes function satisfies 

(4.13) WD x f = 2VU x f + U, 

where U(x) in is a polynomial in x (as are V(x) and W(x)). 

Proof. Starting with the definition (I3.19P we compute 



of 



Buw(u)D x , 

uo<u<u m X — U 



and make the evaluation of the factor 



1 



However we note that 



so that 



1 

x — u (E x x — u)(E x x — u) 
1 1 



1 



x — u (x — Eiu)(x — E u u) 
{x-E+u){x-E-u) n 1 



1 



x — u (EJ x — u)(E x x — u) x — u x — u' 

where we have assumed a symmetric quadratic lattice for simplicity, but point out that our argument could 
be extended to the general case. We deduce then that 

1 



,f = - 



uw(u)H> u - 
l 



{E+u) , B u w(E+u) - 

UQ<U<U<yi % U 



>(E+*uk) 



U)(U ) 



x - E£u<n x- EuUo' 
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after employing the formula for integration by parts. We next use the definition of the D-semi-classical 
weight (|4.ip and perform a sequence of subtractions in the numerator of the integrand 



Ju <u<u m x-EZu W(EZu) x-EZurn x-E u u 

= 2V(x)M x f(x) + 2V(x) { f B(E+u) U ^ E ^ _ Mx f(x) 



2 I D(E+u) 



MEju) V(E+u)W(x) - V(x)W(E+u) , „ 



W{E+u) x - E+u 



Since the third term on the right-hand side of the previous equation is clearly a polynomial in x we shall 
focus next on the second term. The last factor of this latter term can be written as 

M x f(x) - / B(Eiu) = / !*«,(«)* + } 



1 



J Uq<U<U,jx X £J U U J UQ<U<Uyi X Fju 

The last two terms of the above equation are 

> u s ] — — + > [£+ m s -m s ] — { —+Ay[E+u m ) -— Ay(E u u )- 



n s x-E u u s x-EZu s x-EuU<n x-E u u 



E, \ f ^i" 2 "s - . u s~ E u 2 U S 
W{U S ) < T 1 = 
, . , m L x - EZu s x- E u u s 



o<s<m 

where we have used the identity for integrals under the re-parameterisation of the lattice u i— > E~ 2 u and 
introducing the additional boundary terms. Combining these results we have 

1 1 \ K+V- - u s u s - £~ 2 w s 



X 



E£u s x - E u u s 



lAy(£;+ MOT ) 1 ' ^ + ^Ay(E-u ) ^ 



= 1 V x(n s ) [{E+u s E~u s ) ^ E i +E ^ x - 2Us E+ 2 u s -u s u s -E u 2 u s \ ^ 

2 0<T<« ^ " (s-^ttU^^-Satt,) X-Etu s X-EuU s J 

where we have used the ratio identity again. We observe that the summand appears to have two sets of 
simple poles at x — E^u s , however when we compute their residues they vanish identically, implying that 
the summand is in reality a polynomial in x. This leaves us with the two sets of boundary terms to account 
for. We find that their combination is 

_ w(E+ 2 un) [w{x) _ AyiE + Urn)v{x)] + «;(uo) _ Ay{E - UQ)v(x)] _ 

x - EZuyi x - E u u 

The first term vanishes under the upper terminating condition whilst the second is a polynomial in x due to 
the lower terminating condition. In summary this leads us to conclude that W{x)U x f(x) — 2V(x)M. x f(x) is 
a polynomial in x. □ 

We offer the following analog to the concept of a regular semi-classical weight. 

Definition 4.2. A generic or regular D- semi- classical weight has two properties - 
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(i) strict inequalities in the degrees of the spectral data polynomials, i.e. deg x W — M, deg x y = M — 1 
and deg x U — M — 2, and 

(ii) the lattice generated by any zero of (W 2 — Ay 2 V 2 )(x) does not coincide with the lattice generated 
by another zero, i.e. for any x such that (W 2 - Ay 2 V 2 )(x) = then (W 2 - Ay 2 V 2 )(E 2Z x) ^ 0. 

Conversely any weight with data not satisfying these conditions will be termed an irregular D- semi- classical 
weight. Note that the zeros of W are not relevant and no interpretation as a singularity can be placed upon 
a zero of (W 2 - Ay 2 V 2 ){x). 

Proposition 4.2 (|63J). The spectral coefficients W n (x), Q n (x) and 0, n (x) we defined in terms of 'bilinear 
formulae of the polynomials and associated functions 

(4.14) —(2W n (x)-W{x)) = {W + AyV)\p n {y + )q n ^{y^)-p n - X {y + )q n {y^)] 

+ (W - AyV) [p n (y-)q n -i{y+) - p„-i(y-)q„(y+)] , n > 0, 

(4.15) Aye n (x) = (W + AyV)p n (y + )q n (y-)-(W-AyV)p n {y-)q n (y+), n > 0, 

(4.16) ^-(Q n (x) + V(x)) =(W + AyV) \p n (y + )q n ^(y.) + p n ^(y + )q n (y^)] 

+ {W- AyV) [-p n (y-)q n -i(y+) - J>n_i (»_)?„(»+)] , n > 0. 
For the D- semi- classical class of weights these coefficients are polynomials in x. 

Proof. Our proof will be an extension of Laguerre's reasoning for the differential case. However we first need 
to appreciate one important fact. Any polynomial in y_ and y+ is an element of the Ring <C[x] + AyC[x] 
because y± — M x a;± \Ay and M x i is a linear polynomial in x and Ay 2 is a quadratic polynomial in x. The 
starting point will be the definition for q n , but inverted to solve for the Stieltjes function 

a W 1} 

(4.17) / = qn+P "-\ 

Pn 

for arbitrary n > 0. Now consider the formulation of (14. 13)) using this expression 

0=(W- AyV)f(y+) - (W + AyV)f(y-) - AyU, 
= (W - AyV)q n (y + )p n (y^) + AyV)q n (y_)p n ( y+ ) 

p n (y+)pn(y-) 

| (W - AyV)p^_ 1 {y + )p n {y^) - (W + AyV)p^_ 1 {y_)p n {y + ) - AyUp n (y + )p n (y_) 

Pn(y+)pn(y-) 

Now (W - AyV)p^ ) _ 1 (y+)p n (y-) - (W + AyV)p^ ) _ 1 {y-)p n {y+) - AyUp n {y+)p n {y.) G C[x] + AyC[x] and 
furthermore because it is odd under y + <-> y_ this quantity is actually an element of Aj/C[x] only. Conse- 
quently (W - AyV)q n (y+)p n (y-) - (W + AyV)q n (y-)p n (y+) G AyC[x] and we define the polynomial 6„(x) 

by 

(4-18) AyO n (x) = {W + AyV)q n (y-) Pn (y + ) - (W - AyV)q n (y+)p n (y-), 

which is (|4.15[) . To find the other relations we employ the alternative form of & n (x) 

(4.19) A y e n = {W- AyV)p^\{y+)p n {y^) - (W + AyV)p^ ) _ l {y-)p n {y+) - AyU Pn (y + ) Pn (y.), 

and multiply the left-hand sides by the Casoratians 1 — a n [p n -\{y±)p t ^}_ 1 (y±) — p n {y±)p^-i{y±)] so that 
upon re-organising we have 

(4.20) Pn (y+) [(W + AyV)p i ^_ 1 (y-) + &yU Pn (y-) - A y a n Q nP { n 1) _ 2 (y + ) 

= Pn-i(y+) [(w - &yV)Pn(y-) - Aya n e nPn - 1 ( y+ )} , 
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and 

(4.21) p n 0/_) [{W - AyV)p^l 1 (y + ) - AyUp n (y+) + Aya n Q n £ ] _ 2 { V -) 

= Pn-iiV-) [{W + &yV)p n (y+) + Aya n Q n p n ^(y-)] . 
Relation (|4.20l) implies there must exist polynomials 7ri,7T2 € C[x] such that 

(4.22) p n ( y+ )p { r ll 1 (y + )(7r 1 + Ayn 2 ) = p n (y+) [(W + AyV)p^_ 1 { y ^) + AyUp n (y.) - Aya n Q nP ^l 2 (y + ) 

(4-23) = p^lM [(W - AyV)p n (y-) - Aya„e n p n . 1 (y + )] , 

while (|4.2ip implies an identical relation with y + <-> t/_ and Ay i— > —Ay. Therefore we have the pair of 
relations 

(4.24) pW 1 (y + )( 7 r 1 + Ayn 2 ) = (W + AyV)p^_ 1 (y^) + AyUp n (y.) - Aya n Q n p^_ 2 (y + ), 

(4.25) Pn(y+)(7ri + Ayir 2 ) = (W - AyV) Pn (y-) - Aya n Q nPn ^ 1 (y + ), 

and an analogous pair with y + «-» We form the combination of p n ^i{y + ) x (|4.24|) — p^_ 2 (y + ) x (|4.25[) 
and by employing the Casoratian and eliminating the associated polynomials we arrive at the solution 

(4.26) TTi + Ayvr 2 = -a„(W + AyV) Pn -i{y+)q n {y-) + a n (W - AyV) Pn (y-)q n -i(y + ). 
A similar procedure on the other pair yields the partner relation 

(4.27) TTi - Ayn 2 = a n (W + AyV)p n {y+)q n -i{y-) - a n (W - AyV)p n -i(y-)q n (y+). 

The sum and difference of these are simply (|4. 14[) and (14. 16)) with the identifications 7Ti = 2W n — W and 

7T 2 = -0„ — V. □ 

Knowledge of the large x expansions of the spectral coefficients will be important in the ensuing analysis. 

Proposition 4.3 ( 60 ). The spectral coefficients have terminating expansions about x = oo with leading 
order terms 



(4.28) W n {x) = \W+\[W + AyV](^j + \[W-AyV}(^) +0(x M - 1 ), 



n > 0, 



(4.29) e n(a; ) = ^_ [W + AyV] (?±) - -^-[W - AyV] ( V — ) + 0(x M - 3 ), n > 0, 

V-Ay \y-J y+Ay \y^ 

(4.30) n n (x) + V(x) = -^-[W + A y V}(^] - ^-[W- AyV] (— ) +0{x M - 2 ), n>0. 



2Ay L \y-J 2Ay l \y+ , 

Proof. These expansions follow from the substitution of the expansions for the polynomials (|3.34l) and 
associated functions (|3 .35[) into the formulae (|4.14[) - (|4.16[) . One should note that y + /y_ = 0(1) as x —> 
oo. □ 

Corollary 4.2 f[60]). In the generic or regular D- semi- classical case W n ,Q n ,Q n are polynomials in x with 
fixed degrees independent of n, specifically deg x Wn = M , &eg x Sl n = M — 1 and deg^O™ = M — 2. 

Proof. This follows from the previous proposition after recalling that y± = 0(x) as x — > oo. □ 

The expressions for the spectral coefficients (|4. 14114. 161 can be inverted yielding a system of coupled first 
order linear divided-difference equations for the polynomials and associated functions with respect to x. 

Proposition 4.4 ( 60, 63J). The polynomials with a J}- semi- classical weight satisfy 

(4.31) W n D e„M lPn _i, n>0, 

for certain spectral coefficients, W n (x) 1 Q n (x),fl n (x). The partner equation to 31\ ) is 

(4.32) W n B x p„_i = -(n„ + 2V)M xPn _ 1 + a n Q n -iM xPn , n>0. 
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The associated functions satisfy divided- difference equations 

(4.33) W n B x q n = (ft„ + 2V)M x q n - a„G n M x q n - U n > 0, 

(4.34) W„D e n -iM x q n , n>0, 

identical in structure to ^4.31\ ) and ^4.3<fy respectively. The coupled system fr4.31\4-32\4-33\4-34\ ! can be 

written in matrix form as the spectral divided- difference equation 

(4.35) D^M^A^^).^^^ _ Q ^f% x) )M x Y n ( X ), n > 0, 
with A n termed the spectral matrix. 

Proof. Our derivation of the spectral divided-difference equations will employ a simple method - that of 
inversion of the relations for the spectral coefficients (|4. 14114. f 6]) . Consider the following equivalent forms of 
those relations 

(4.36) 2W n - W + Ay(Q n + V) = a n [(W + AyV)p n (y+)q n - 1 (y-) - (W - AyV)p n - 1 (y-)q n (y + )} , 

(4.37) 2W n - W - Ay(n n + V) = a n [-(W + AyV)p n ^ 1 (y+)q n (y^) + (W - AyV)p n {y-)q n - l {y + )\ , 
(4-38) AyQ n = (W + AyV)p n (y+)q n (y-) - (W - AyV)p n (y-)q n (y+). 

For our first task, that of deriving (|4.31|) . we construct the combination p n (y+) X (|4.37j) — Pn(y-) x f|4.36|) + 
(Pn-i(y+)+Pn-i(y-))a n x (|4.38[) . After effecting some initial cancellation we can then employ the Casoratian 
relation (|3.26[) to make further simplification and our final result is precisely (|4.31[) in finite difference form. 
The approach taken for the other relations (I4.32[) - (I4.34|) is the same. For (|4.32l) we use the combination 
Pn-i(y+) X^M-Pn-i(y-) X$M-(Pn(y+)+Pn(y-))a n x (|3381 1 UI1U r while for we use q n {y+) x 

(14371) - q n (y_ ) X (OBI + (g»_ x (t/ + ) +g w -i (y-))a n X flSlJl and for Sttity we use g n _i(j/ + )x ^M-Qn-xiV-) X 
(1132) - (Qn(y+) + q n (y-))a n X $m\n^n-V D 

The initial values of the spectral coefficients are specified by the weight data W, V, U and the initial norm 

7o 

(4.39) W-i=0, W = W, e_ x =0, 9o = -7o^ ^o = 0, Ox = -2V - ^{M x x - b )U. 

Remark 4.4. An alternative line of reasoning to achieve the results above, and in fact to generalise them 
beyond the ©-semi-classical class, is the following. By positing an explicit quadrature formula for the non- 
linear lattice, one upon which the orthogonality condition is founded, one can derive (|4.31[) , (|4.32p and (|4.35j) 
assuming the existence of certain moments of "log-derivative" of the weight H) x w/M. x w. A byproduct of this 
approach is that one also constructs an explicit quadrature formula for the spectral coefficients, albeit one 
which involves the polynomials themselves. 

The spectral structure in Proposition ^. 4l has been derived before for OPS with respect to general classes of 
weights on the linear lattice by [29l EH SU [68] and on the q-linear lattice bv [55l 167] |4S ] 154 ] [20 ] l32l l46l 170] [79] . 
This has also been done for specific weights on the q-linear lattice in [72] [15] and for OPS on the unit 
circle subject to (/-difference relations in [35] [T3J. In addition to the studies by Magnus the works |2"5] 
have investigated the spectral structures for general SNUL and arbitrary deg x W,deg x V . Characterisation 
theorems for the classical OPS of the Askey table have been made by [22] [T| and others using an equivalent 
formulation to the spectral structures on general SNUL, but of course restricted to the deg^VF = 2, deg x V = 1 
case. 

Compatibility relation between the matrix recurrence relation and the spectral divided-difference equation 
imposes the following conditions. The Cayley transform (l — ^Ay A) (l + \ Ay A) will figure prominently 
is the ensuing analysis. 

Proposition 4.5. The spectral matrix and the recurrence matrix satisfy 

(4.40) A n+1 -M x K n -M x K n - A n =H x K n - \Ay 2 A n+1 -H x K n ■ A n , n>0, 
or equivalently 

(4.41) K n {y+) (1 - \AyAnY 1 (l + \Ay A n ) = (l - \Ay A n+1 ) _1 (l + \Ay A n+1 ) K n (y_), n > 0. 
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Proof. The finite difference form of (|4.35p is 

(4.42) (1 - \AyA n ) Y n (y + ) = (l + ±AyA n ) Y n (y_), 

which demonstrates why the Cayley transform arises. The second relation (|4.4ip is found by computing 
Y n +i(y+) in two different ways from Y n (y-) corresponding to the different orders of the operations n4 n+1 
and y_ ^ y+. The first relation (|4.40l) is derived by computing H x Y n +\ in terms of M. x Y n in two different 
orders. □ 

Compatibility of the three term recurrence relation and the spectral divided-difference equation equation 
implies the following result. 

Proposition 4.6 ([60, 63 ). The spectral coefficients arising in Proposition \4-4\ satisfy recurrence relations 
in n, 

(4.43) W n+1 = W n + \Ay 2 Q n , n > 0, 

(4.44) 0„ +1 + fi n + IV = (M x x - b n )Q n , n > 0, 
(4.45) 

(w n n n+1 ~ w n+1 n n )(M x x - b n ) = -\A y 2 n n+1 n n + W n W n+1 + a 2 l+1 w n e n+1 - a 2 n W n+1 Q n -i, n > 0. 

Proof. We offer a proof in the spirit of our previous proofs. For the first relation (|4.43l) we construct 
2W n +i — 2W„ — ^Ay 2 Q n and substitute the expressions (|4.14p and (14. 15)) for W n and 8„ respectively. 
After collecting terms proportional to W ± AyV we observe that the coefficients of these terms separate 
into two parts, in which the three term recurrence relation can be employed. The result we find is that 
they identically vanish. The second relation (|4.44|) can be found in an identical manner starting with 
Ay [Q n+ i + fi„ + 2V - {M x x-b n )Q n \. □ 



The relations of Proposition 14.61 constitute the analogs of the Laguerre-Freud equations and have been 
studied in the special case of linear lattices by [29] [65] 24] [68] and in the g-linear case by [55] [54] [TO, 32 . For 
the situation of OPS on the unit circle with q-difference operators some of these relations can be found in 
[48] . Following on from the Magnus studies the works [9] [28] have investigated these compatibility conditions 
for general SNUL. Together with Proposition 14.41 this result allows for the construction of ladder, or raising 
and lowering operators and these have been found in the q- linear lattice by [50] SS] . 

A consequence of Proposition 14.61 is the following. 

Proposition 4.7 ( 60 ). The spectral coefficients satisfy the bilinear recurrence relation 

(4.46) W n (W n -W) = \Ay 2 [n n (Q n + 2V) - a 2 n Q n ^e n ] = -\Ay 2 dct (Jq^ Jn^v) ' n ~ °' 

Proof. Two methods of proof are available here. The first method starts from the observation that the above 
bilinear relation can be expressed as 

(4.47) [2W n -W + Ay(fl n + V)] [2W n -W- Ay{Q n + V)] + a ? 2 l A 2 / 2 e n 6 I1 _i = W 2 - Ay 2 V 2 . 

Now we have evaluations of the factors appearing in the left-hand side of the preceeding equation from the 
proof of Proposition 14.21 namely (|4.261 14.271 14. 18[) . in terms of products of a polynomial and associated 
function. Substituting these evaluations into the left-hand side, effecting considerable cancellation of terms 
and employing the Casoratian relation p.26p the final result simplifies to the right-hand side. 

The second method uses the compatibility relations (14.431 14.45P and constructs an integral for this coupled 
system. First one multiplies the left-hand side of (|4.44l) by the left-hand side of (14.45)) and cancels out the 
common factor of Ma; — b n in the resulting equation. By judicious use of (|4.43|) one can effect the cancellation 
of four terms in this equation and a recasting of the terms linear in Q n . Consequently we arrive at the 
following equation 



(4.48) V 1 ~ V n J r ' =— 2 {W n+ ,-W n ) + V 2 



1 1 



w n+1 w n 



2 ®n+l®n 2 ®n®n- 



W n+1 " W n 



(n n+1 + v) 2 (n n + vy 

W n+1 W n Ay 2 

which is manifestly a perfect difference equation in n. By summing this and employing the initial value 
evaluations (|4.39j) we have (|4.46|) . □ 
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The above Proposition must be augmented with the following matrix identities. 
Corollary 4.3. The matrix factor appearing in Proposition ^. 5\ has the determinant evaluation 

W =F Aw V 

(4.49) det(l ± % Ay An) = T y - , n > 0, 

w n 

and its inverse is 

(4.50) (l±±AyA n )-i= \ (W n ^Ay(V + 2V) ±^ y a 9 \ > 
V y V 2 y ; WTAyV\ T^Aya n Q„-i W n ±\AyVL n ) 1 

Proof. From the form for the spectral matrix given by (I4.35|) we have 

det(l ± \Ay An) = [(W n ± |Ayfi n )(W n T |Ay(Q„ + 2V)) + \Ay 2 a 2 n Q n Q n ^] 

= [W 2 T AyFW„ - ±Ay 2 [fi n (n„ + 2F) - a 2 i e n e n _ 1 ]] 

_ Wt AyV 

W n ' 

where we have employed the previous Corollary in the last step. Constructing the matrix inverse given on 
the left-hand side of (|4.50|) and using our determinant formula we find (|4.50[) . □ 

Consequent to the results of Corollary 14.31 we require the matrix product found in (|4.41l) . 
Corollary 4.4. The matrix product appearing in ^4.41\ l has the evaluation 



n > 0, 



(4.51) (l-±AyA n ) l (l + ±AyA n ) 

= 1 (2W n -W + Ay (a n + V) -Ay a„9„ 

W + Ay V \ Ay a„6„_i 2W n - W - Ay (fi n + V) / 

= - • <«-'<»->■ -*<»-)) - ( *^,) « (p-fe-), -*&-» ■ 

Proof. Using (|4.50j) we multiply out the matrix product and employ the bilinear identity (|4.46j) to simplify 
the diagonal elements, resulting in (|4.51j) . □ 

This result motivates the following definitions 

(4.52) 2»± :=2W n -W± Ay(Q n + V), n > 1, 

(4.53) T+ := Aya n Q n , T_ := Aya n 6 n _i, n > 1, 

whilst for n = we have W±(n = 0) := W ± AyV, T+(n = 0) := -Aya ^U, T_(n = 0) := together 
with 



(4.54) A* n := 



The variables S2JJ± are essentially the variables U n ,X„ employed by Magnus in 63 respectively whilst the 
variables %± are directly related to his Y n , Z n variables. In these new variables the bi- linear relation (|4.46[) 
becomes 

(4.55) det A* n = 2U+2U_ + T+T_ = W 2 - Ay 2 V 2 , 

which will be the more useful form. 

It is of advantage to write out recurrence-spectral compatibility equations in terms of 2U±, n , T± )7l , where 
we append a subscript to indicate the dependence on the index n. 
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Corollary 4.5. Solving ^4-41$ for A* l+1 we deduce 

(4.56) a„o„+iT +in+ i = -a n (y + - 6„)2U +i „ + a„(y_ - 6 n )2U_ >n + (y+ - 6 n )(y- - & n )^+,n + 

(4.57) a„T 

(4.58) a„2U +: „ + i = a„2B_,„ + (y+ - &„)1+,„, 

(4.59) a„2U_, n+ i = a n 2» + ,„ - (y_ - b n )1+, n . 

A significant result follows from the observation that the right-hand side of (|4.55p is a polynomial in x 
with fixed degree independent of n and in fact contains no dependence on n. If we denote one of the zeros 
of the spectral polynomial W 2 — Ay 2 V 2 as Xj then we can apply the equality for n H> n + 1 and use (|4.57|) 
to deduce 



(4.60) 



This then allows us to draw the following conclusion. 

Proposition 4.8. The variable X+ ! , l +i(x :) ), when written in terms of the variables at n, factorises in the 

following way 

(4.61) 

- a n a n+ i1 +tn+ i(xj)% +:Tl (xj) = [a n W- in (xj) + (y +ij - b n )% +in (xj)] [a n W +in (xj) - - b n )% +!n (xj)] . 

The coupled pair of first order divided-difference equations (|4.3ip and (I4.32[) imply a second order equation 
for one of the components, say p n . 

Proposition 4.9 ([60, 8 ). The D- semi- classical orthogonal polynomials or associated functions satisfy a 
second order divided- difference equation in the following two equivalent forms - either 



(4.62) 



(W_n 



7^ + * VD a 



e„ 



Wn 
On 



w-w n 
A y 2 e n 



( W- w n 

a 
e 



or alternatively 

fW + AyV\ 



(4.63) £+ 



-E. 



EZ 



_ f W-AyV \ 
x \ AyQ n ) 
2W n -W + Ay(Sl n + V) 



AyQ n 



(KVPn 

EZ 



2W n -W- Ay(Q n + V) 
AyQ n 



0. 



Proof. Starting with (14.31)) and (|4.32p . we can employ (|4.46l) to write an equation for each of M x p n _i and 
J3 x p n -i in terms of linear combinations of M. x p n and H x p n . Then utilising the fact that M. x H) x p n -i = 
DiMiP„-i we arrive at (|4.62[) . The equation (|4.63[) is the nodal equivalent of the former. □ 

This concludes our discussion of the spectral structures for a general quadratic lattice and we now turn 
to the simplest explicit example on the q-quadratic lattice. 



5. M = 2, L — CASE AND THE ASKEY-WlLSON POLYNOMIALS 

Here we will employ the theory of Section [4] and explicitly compute the spectral coefficients for the Askey- 
Wilson system itself because it serves both as an instructive example for the theory of the previous section 
and clarifies some confusion present in the literature. Virtually all of our results presented here have been 
found by earlier studies - firstly by Askey and Wilson [3] , then most notably by the Soviet school of Nikiforov, 
Suslov and Uvarov [761 1731 [771 1741 175) who were primarily concerned with hypergeometric type OPS on non- 
uniform lattices, and in the 1988 work of Magnus [60 whose work and intent is most similar to the spirit of 
our own. 
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(5.1) 



We recall the Askey- Wilson weight [J] itself has degrees 2N = 4 or M — 2 and L = with 

(z ±2 ;q)oo 



){x) = w(x; {ai, . . . ,a 4 }) = — 



Bin(8)U.j=i( a 3 z±1 ><l)<: 



-, G = (-l,l). 



Let the j-th elementary symmetric polynomial of ai, . . . , 04 be denoted by Oj for j = 0, . . . , 4. 

Proposition 5.1. Assume that g ^ 1 and | <7 1 / 2 cj ^ | 7^ 1 /or j = 1, ... ,4. XTie spectral coefficients for the 
Askey-Wilson OP system are 



x + 1 + q 1 <7 2 + q 2 0"4, 



(5.2) W n (x) = (1 + <T n )(l + ff 4 g ,l - 2 )(x 2 - 1) 
(5.3) 

-Wax + (-2 + q- n )q-V 2 a 3 + (-2 + q n )q-^l 2 o x o A + q n - 7 ^a 3 a 4 



n n (x) + v(x) = 2 qn ,y 4 *" x-< q ~ 



g l/2 _ g -l/2 



( ? l/2 _ g-l/2)( g -n _ g n-2 CT4 ) 



(5.4) 



_ g ~-3/y 4 _ g -n-l/2 
^„^Xj — "t ^ 2 _ 1/2 i 



valid for n > 0. TTie three-term recurrence coefficients are found to be given by the standard expressions 



(5.5) 
and 



1 (l-a")(l-a4g"- 2 )n fc > i (l-« J «fcg 
4 (1 - a 4 q 2n - 3 ){l - a 4 q 2n - 2 ) 2 {l ~ ^g 2 ™" 1 ) ' 



(5.6) 6„ = [a x (g + a A (q 2n - <f - g"" 1 )) + a 3 (l - g" - + ^g 2 "" 1 )] 
for n > where we assume oq = 0. 

Proof. From the weight (15 . 1 [) we compute the spectral data 



„n-l 



2(l-(74g 2 ")(l-CT4g 2 ' l - 2 ) : 



(5.7) 



W ± AyV = z^ 2 JJ(1 — a,^ 1/2 z ±1 ), 



from which we deduce 

(5.8) J¥(a;) = 2(1 + < 7 4 <r 2 ) ; r 2 



a; — 1 + q 1 a 2 - q 2 o"4, 



(5.9) 



V(x) = 2 



q V4 — 1 g 1 / 2 (Ti — g 3 / 2 o"3 



,1/2 _ -l/2' 



g 



1/2 _ fl -l/2 



We will see that W 2 — Ay 2 ^ 2 will play a significant role and therefore define another set of elementary 
symmetric polynomials by 

(5.10) W 2 - Ay 2 V 2 = K 2 [a; 4 - e lX 3 + e 2 x 2 - e 3 x + e 4 ] . 

We note the evaluation K = \q^ = 4q~ 1 v / (74. We parameterise the spectral coefficients in the following way 

(5.11) 2W n - W = w 2 x 2 + wix + w , 

(5.12) e„ = o7 + , 9„_i=n7_, 

(5.13) fl n + V = V\x + vq, 

where we know from (|4. 28114. 3"0|) that the leading coefficients are 

K , „ 4 .„ , K q n +»-q- n -» 



(5.14) 



„n+ju+l/2 _ -n-/j-l/2 

"^t^" + *—"). "' = i»- -+ = ^ g i/2_ 9 g -i/2 — • 
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From the fundamental bi-linear relation (|4.46l) we get a system of quadratic polynomial equalities 

(5.15) w\ - Av\ = K 2 , 

(5.16) 2w 1 w 2 -2Av v 1 =-K 2 e 1 , 

(5.17) 2w w 2 +wj- A(v 2 - v{) + a 2 n Aw + xu- = K 2 e 2 , 

(5.18) 2w wi + 2Av vi = -K 2 e 3 , 

(5.19) wl + Av 2 - a 2 n Aw + W- = K 2 e i7 

where A := (q 1 / 2 - q- 1 / 2 ) 2 . Now (15.1515. 1715. 19[) imply wl + (w Q + w 2 ) 2 = K 2 (l + e 2 +e 4 ) while (|5.18I5.16|) 
imply 2wi(wq + w 2 ) ~ — K 2 (e± + e 3 ) . Forming the sum and difference of these two later relations we conclude 
that 

, eort , W{l)-W(-l) W(l) + W(-1) 

(5.20) wi = 2 — = -w 2 + 2 — -■ 

Using (|5.16[) for example, along with the above solutions, we find 

(5.21) ^e 1+ ^(l)-^(-l)) 
v ' 2Av x 

Therefore we have succeeded in relating the sub-leading coefficients in terms of explicitly known quantities 
and after simplification we arrive at (|5. 2115. 31) . Finally, using (|5.19[) , we can solve for a 2 n and after observing 
that e4 can be expressed as 



q-» W(l) + W{-1) 



2 



1 



qf + q-f W(l) -W(-l) 
ei H 



K 



2 



we can factorise the resulting four terms and arrive at (|5.5[) . To find (|5.6|) we start with (|4.44l) and use our 
previous results for (|5 ,4[) and (|5 .3|) . We note the terms linear in x cancel identically, as they must, and after 
some factorisation we deduce (15.61) . □ 



It is the coefficients (|5.2II5.3|) along with first order divided-difference equation ()4.31|) that constitutes the 
structural relation for the Askey- Wilson polynomials, that is the analog of first order difference or differential 
relations for the classical orthogonal polynomials. The divided-difference relations reported in [581 are all of 
second order. 

It is an easy task to evaluate the Askey- Wilson integral [H [56] as a g-factorial and therefore compute the 
moments. The Askey- Wilson integral is defined by 

5-23 h{a\, a 2 ,a 3 ,a 4 := / ^—7^3 — : t: — — . 

y T 2^zn 4 =1 (a J z ±1 ;q) 0O 

with \a,j\ < 1 for j = 1, ... ,4. 

Our method is to apply the general system of moment recurrences to the case M = 2 and we find that it 
coincides with the recurrence of Kalnins and Miller |53) and Koelink and Koornwinder [57) . 

Proposition 5.2 ( 97J). The Askey-Wilson integral satisfies the two-term linear recurrence 

(5.24) (cr 4 - 1)12(901,02,03,04) = (aia 2 - l)(aia 3 - l)(aia 4 - 1)12(01,02,03,04), 
which is solved by 

(5.25) /2(ai,02,a 3 ,a 4 ) = 2- — ^ 4 ' q ^ 00 

(<?; 9)oo H k>J \aja k ; q)^ 

Consequently the moments are given by 

t e og n / \ (aia2,aia 3 ,aia4;g)„ 

(5.26) m ,„(ai) = n -. r 12(01,02,03,04) 

(0-4; q)n 

(5.27) =2n (oV 4 ;g)oo 



(q n aia 2 ,q n aia 3l o n aia 4 , a 2 a 3 , a 2 a 4 , 0304, q; q) c 
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Proof. We only require the k = case of (|4.11j) with a — a\ and from (|5.7p we compute the relevant 
coefficients as 62,1(0,) = q~ 1 a~ 1 (l — 04) and <52,o(a) = <Z _1 ( a + a -1 )^ — 1) + q _1 (o"i — 03). The solution of 
the recurrence (|5.24l) follows from the arguments given in |S7]- □ 



The specialisation of the second-order divided-difference equation (|4.63p with the M = 2 spectral coeffi- 
cients of Proposition 15 . 1 1 is given by 
(5.28) 

Ogl^Z^fj [ {E +fp n - Pn ]+U^Z^l[ {E -f pn _ Pn ] +q -n-l {1 _ q n ){ n ai _ ){z 2 _ 1)pn = Q 
qz z — 1 z z — q 

This is soluble in terms of basic hypergeometric functions and their polynomial solutions have been given by 
[6"U] or using factorisation methods in [TJHHHO] as Askey- Wilson polynomials. The Askey- Wilson polynomials 
have an explicit form as a balanced 4^3 function [56 with the manifest symmetry under z O z^ 1 and are 
given in the monic form by 

(K on x / \ [a 1 a 2 ,a 1 a Sl a 1 a i ]q) n ( q~ n a^q 1 ^ 1 , aiZ, a^ 1 

(5.29) TT n (x) — —- — — = r — 4093 ;q,q 

(2ai)"(g"- 1 cr 4 ;g), i V aia 2 , Oia 3) a lttA 



or the alternative form |34| which is manifestly symmetric under permutations of 01,02,03,04 

fK on\ t \ to ■,-n( a i z ^ a 2Z,a 3 z 7 a 4 z 7 a 4 q- 1 ;q) n _ 2 _ n _ x _ a _ x _ x q 2 ~" 

(5.30) tt„(x) = (2z) p-= — — -: — sW 7 (q z ;q ,a\Z ,a 2 z ,a 3 z ,a 4 z ; ). 

(z z ;q)„(cr 4 q , q)2n C4 

From Eq. (2.6) of [82] or equivalently from Eq. (4.18) of 45 with a = and Eq. (111.23) of [31] we 
deduce the expression for the Stieltjes function as a very-well-poised sWj 
(5.31) 

tt \ M« 1 o-A\q)oo {1-qz 2 ) , -2 1 -1 1 -1 Q -1 1 -1 -1 \ 

/W = 7 s — rf 7 \ rf ^ rTTsWMqz 5 — z ,—z ,—z ,—z ,q;q,q a a), 

(95 < ?)ooll fe>J (ajO fc ;(j) 00 zH 1 < J <4 1 1 ~ a J z J °i a 2 a 3 a4 

which exhibits the parameter symmetry. This expression is the one valid on the second Riemann sheet of 
the cut s-plane or exterior to the unit circle in z, \z\ > 1. There are evaluations of q n (x) in [82] and [45] 
that are simple generalisations of the above formulae but we do not need to discuss them here. 

Having an explicit form for the Stieltjes function we have two tasks at hand - to verify that is satisfies 
the inhomogeneous divided-difference equation (|4.13p and the large x expansion formula (|3.46[) . 

Proposition 5.3. The Stieltjes function given by Eg. \5.3f]) satisfies the inhomogeneous divided difference 
equation |^._?ff| ) with the constant 



(5.32) U 
Proof. We require an alternative form 



_8_7r (q 1 cr i \q), 00 

a-l {q\q)oo\~[ k>:j {aja k ]q) 



q (74 q 2 oiz 

;q)c 



(5 33) f(x)—4-7r qai ' ' a 2 z 7 a 3 z 7 oaz' qa x z' a 4 



04 f g ; 1 8 1 .g\ ) (0203,0204,0304,012: 1 ,a 2 z 1 ,a 3 z 1 ,a 4 z 1 ,a 1 z;q) oc 
a 2 a 3 ' 0204 ' a 3 a 4 ' 

, ai ai ai q 2 
(q cr 4 ,q— ,q— ,q— , — ;q)oo 2 1 

, 020304^4 „, 1 2 -1 q \ 1 

+ 7 2 — \ — 1 \ — T\ T71 —r:&W 7 (a l ;a 1 z, a\z , a x a 2 , aia 3 , aia 4 ; — ) > , 

(q, gafjgjoo ll k>j {aja k ; q)oo(l - aiz)(l - a\Z L ) a A j 

which makes the z z~ l symmetry manifest in the second term only. This is derived from the previous 
expression by utilising Eq.(III.37) combined with Eq.(III.23) of [31], and Eq.(III.36) for the term which 
specialises. We are now in a position to verify that (|5.33[) satisfies (|4.13p . This relies on the fact that (|4.13|) 
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with the data (|5.7[) is precisely the contiguous relation, Eq. (2.2) of [45] . which in our context states 



t roA\ 2^=2, 3,4( 1 1 a 3 Z ) U r , 2 1/2 -1/2 -1 1 \ 

(5.34) z -j- 8 W 7 {a 1 ;q ' diz, g 1 aiz , aia 2 , aia 3 , aia 4 ; — ) 

1 — (7 i / z aiz (74 

-2111=2,3,4(1 ^ 9 ^dj-z) . . 

- z r-pr : ? ,W 7 {a 1 ]q 1 a x z , q 1 a x z , aia 2 , aia 3 , aiar, — ) 

1 — q L ' z aiz 1 (74 

= g~ 3 / 2 a 2 a 3 a 4 (l — )(z — z~ 1 ) & W 7 {a\;q 1 l 2 a,iz, q 1 l 2 a\z~ 1 ,a\a 2 , aia 3 , aia 4 ; — ). 

(74 (74 

We have also used the specialisation formula when aq = be 

Tjr , , , , , 9 2 a 2 n (ag,ag/de,ag/d/,ag/e/;g)oo 

(5.35) 8^7(0; b,c,d,e,f;q, ) - 



bede/ (ag/ d, ag/ e, ag/ /, aq/ de/; 9)00 ' 

which applies to the right-hand side of the previous equation. A consequence one can draw from this 
calculation is the explicit evaluation of the constant spectral coefficient. □ 

Remark 5.1. The constant U can be found in another way. This is by noting the initial condition for $0 
is given by (|4.39[) and that is given purely by the leading order term in (15. 4[) . so that knowing the 
normalisation 7q allows us to make the evaluation. 

Proposition 5.4. The Stieltjes function possesses the explicit large x moment generating function formula 
rr oa\ ti \ 47rai((74;g)oo / g,aia 2 ,aia 3 ,aia 4 

(5.36 fix) = — 7 r — -. r — TT4V3 -1 ; 9) 3 

(<?; ?)oo rifc>j( a i a fe; 9)oo(l -a 1 z)(l-a 1 z- 1 ) V oaiz,qaiz i ,cr 4 

A7t(a 2 a 3 a i z^ 1 ,qz^ 2 ;q) 00 1 / a 2 z _1 , a 3 z _1 , a 4 z _1 



3^2 „ „ „ _i _ 2 ;<?,<? 



z(a 2 a 3 ,a 2 a4,a 3 a 4 ,a 2 z \a 3 z x ,a±z 1 \q)oo <t>oo{x\ai) \ a 2 a 3 a 4 z ,qz~ 

Proof. By applying the transformation formula (III. 36) of 31 to the gWy function in (|5 . 33|) we get a sum 
of two 4<p 3 functions one of which reduces to a 3 tp 2 function. After some simplification we get (|5.36[) . We 
observe that the 4 ip 3 term is precisely the second term of (|3.46j) given our formula for the moments (|5.27l) . 
This allows us to conclude that 

, . , . 47r(a 2 a 3 a 4 z" 1 ,gz" 2 ;q) 0O / a 2 z _1 , a 3 z _1 , a^z' 1 

5.37 /oo [x )= -; — — _ 3 y 2 ! 2 ;<?,<? 

z(a 2 a 3 , a 2 a 4 , a 3 a4, a 2 z 1 ,a 3 z 1 ,a 4 z , g)oo \ a 2 a 3 a 4 z , gz 



□ 



6. Deformation Differences 



In this section we discuss the deformation structures based upon a nonlinear lattice in the u variable (or any 
number of variables for that matter) and the corresponding divided-difference operators. We emphasise that 
many of the results for the spectral structure have a parallel result in the deformation structure although 
there will be crucial differences in the details. Denote the deformation variable u and its forward and 
backward shifts by v± — E^u and their difference by Aw := v + — u_. The reader should be aware that the 
deformation lattice does not have to be the same type as the spectral lattice. Our D-semi-classical weight w, 
defined by Definition 14. 1[ acquires an additional dependence on a deformation variable u and furthermore 
satisfies a deformation divided-difference equation with respect to u. 

Definition 6.1. The deformed D- semi- classical weight satisfies 
(6.1) RB u w = 2SM u w, 

or equivalently 

(6-2) ^ = £±1^). 

w(x; u_) R — AvS 

for polynomials S{x; u), R(x; u) irreducible in x and u £ J. In addition we assume R ± AvS =^ for all 
x e G,u e J and that the deformed OPS exists, i.e. that r y n (v±) 7^ for all n g Z>o- We also require the 
condition j (v+)(R - AvS) + j (v-)(R + AvS) ^ 0. 
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There is a very simple motivation for this relation. In the process of extending the weights beyond those 
of the classical cases in the Askey table by increasing the degrees deg x W > 3, deg x V > 2 the weights 
acquire additional parameters which if suitably chosen sit in a completely reflexive or symmetric way to the 
spectral variable and therefore will satisfy an analogous linear, homogeneous and first order divided-difference 
equation in any of these new parameters. 

There are conditions imposed on the weight given that it satisfies both (|4.I[) or (|4.2[) and (|6.1[) or (|6.2[) . 

Proposition 6.1. The spectral data polynomials W(x;u),V(x;u) and the deformation data polynomials 
R(x;u), S(x;u) satisfy the compatibility relations 



(6.3) 



25 

C ~R~ 

or alternatively 



2V 



1 - ^Ay 2 Av 2 nJ—H u % 



R W 



25 



R 



2V 



W 



2V 



\Ay 2 D x — M„— — jAv 2 B u — M x — 



W 



25 



R 



(6.4) 



W + AyV, ,R + AvS, , W + AyV . ,R + AvS, 

— — {y-;u) = — 7-77(2;; v-)— -—{y+;u). 



R-AvS' 



W-AyV y 1 ^'R-AvS^ ' ' W — AyV 

Proof. From the fact that the weight satisfies the over-determined system of first order divided-difference 
equations (|4.1[) and (|6.1|) we first compute that 



(6.5) 
(6.6) 



1 - ^Ay 2 Au 2 I 
I - jgAy 2 Av 2 l 



using (|2.1ip . Now we compute 
(6.7) 
(6.8) 
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x — 
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2* 
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D u l 


) x w = 




\w — 



2V_ 

25 



±Av 2 l 

4 



jAy 2 B x 



— P 


. 2S 
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^R 


— IV 
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u W 



i x m u w, 



,™1 

' W 

25, 



2V 
, — I 
W 

25^ 



» x m u w, 



R R 

using (|2.I0p and use the results of the previous set of equations to derive two independent relations linking 
J}JS) x w and MnMjW. Comparison of these two latter relations leads to (|6 . 3|) . 



To establish (16.41) we note 



W + AyV 

(6.9) w(y+;v + ) = — — -(x; v+)w(y_; v+) 



whereas 
(6.10) 



w(y+;v + 



W -AyV 

_ R + AvS 
' ~ R-AvS 



W + AyV, ^ R + AvS , 



(y + ;u)w(y+;vJ 



W -AyV 

R + AvS 
~ R-AvS 



R-AvS 



W + AyV , 



W - AyV 



□ 



Proposition 6.2. As a consequence of Definition ROl the Stieltjes transform satisfies the inhomogeneous 
equation 



(6.11) 



RB u f = 2SM u f + T. 



The deformed D- semi- classical class of orthogonal polynomial systems are characterised by the property that 
R(x;u), S(x;u) and T(x;u) in f6.11)) are polynomials inx. 



Proof. From the definition (|3.19[) and (|6.1[) we compute 



(6.12) 



D u f(x;u) 



l u w{y;u) 2S{y\u) 



x-y R{y\ u) ' 

Now we observe that the rational function 25(y; u)/[{x — y)R(y; u)] has the partial fraction expansion 

2S{y;u) _ 2S(x;u) | 2S( Xj ;u) 1 

(x-y)R(y;u) (x-y)R{x;u) ' B! {x 3 \v){x - x s ) y - x{ 
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Consequently 
(6.14) ] 



(6.15) 



Kf(x;u) 



2S(x;u) 
R(x; u) 

2S(x;u) 
R(x; u) 



■j.w(y;u) 



x-y 



E 



2S{xf, u) 



R'(xj;u)(x- xj) 



E 



R'{xj] u)(x — Xj) 



l u w(y; u) 
y-Xj 



j,iJ(xj-)=0 

We conclude that 1 1[) follows with deg x T < deg x R — 1. 

Proposition 6.3. Compatibility of H6.11\) and l{4-13\ l implies the following identity on U and T . 
~{W + AyV){x;v^ 



□ 



(6.16) Ay 



(W + AyV)(x; u_ 
= At; 



-(R + AvS)(y_;u)U(x; vJ) -{R- AvS){y-,u)U(x; v + ) 



(W + AyV)(x; v + )T(y.;u) - (W - AyV)(x; v + ) { * ^^1^- T(y + ; u) 



(R-AvS)(y+;u 
Proof. We begin by defining the ratio 

(W + AyV)(x;v+) (R + AvS)(y-;u) _ (W - AyV)(x; v+) (R - AvS)(y^ ; u) 
{ ' ' X (W + AyV) (x;v-) (R + AvS) (y+ ; u) {W - AyV)(x;v-) (R - AvS)(y+; u) ' 
by virtue of (|6.4I) . Consider f|4. 1 3[) in the form 

(6.18) (W - AyV){x- u)f(y+;u) — (W + AyV)(x; u)/(j/_; u) - AyU(x; u) = 0, 

and construct the combination (W + AyV)(x; v+)/(W + AyV)(x; v-)(R + AvS)(y-: u) x lj6TT5]) fa i-> t>_) - 
(R — AvS)(y-;u) x (|6.18[) (u i— > v + ). Using the above identity we find the two terms containing / possess 
the factors {R + AvS)(y_ ; «)/(y_ ;v-)-(R- AvS) (y_ ; u)/fo_ ; v+) and (i2 + AvS) {y+ ; u)/(j/ + ;«_)-(«- 
AvS)(y + ; u)f(y + ; v + ), into which we can apply (|6.1ip in the form 

(6.19) (R - AvS){x; u)f{x; v + ) -(R+ AvS){x; u)f(x; u_) - AvT(x; u) = 0. 

Then (|6.16p immediately follows. □ 



We can extend the notion of a generic or regular D-semi-classical weight, given in Definition 14.21 to the 
deformed situation by the following definition 

Definition 6.2. A regular, deformed D- semi- classical weight is one that satisfies the requirements of the 
Definition 14.21 and deg^i? = deg x S — L and deg^T = deg^i? — 1. Clearly M in Definition 14.21 and L are 
related, depending on the specific case on hand. 

The analog of Proposition 14. 21 is the following. 

Proposition 6.4. Let the deformation coefficients R n (x;u),T n (x]u),'E. n (x;u),^ n (x;u)^ n (x;u) be defined 
in terms of bilinear formulae involving products of the polynomials and associated functions by 



3.20) — Rn = (R + AvS) 

tin 



1 



a n (vJ) 



+ Pn-i(; v+)q n (; vJ) - p n {; w+)g„_i(; «_) 



{R - AvS) 



.21) 



Av 



T n = (R + AvS) 



a n {v+) 

Pn-i{;v + )q n (;v^) -p n {;v+)q n -i{;v-) 



+ Pn-i{;v-)q n {;v + ) -p n (;v_)q n _x(;v + ) 



(R - AvS) 



1 



(6.22) 



Av 



(R + AvS) 



a n (v-) 



a n {v + ) 

+ p n - 1 (;v + )q n (;v-) + p n (; v + )q n -x(; u_) 



+ p n - 1 (;v^)q n (;v+) + p n (; V-)q n -i(; v + ) 



(R - AvS) 



a n (v + ) 



Pn-i(iv-)q n (;v + ) -p n (;v^)q n _i(;v + ) 



2S 



(6.23) 

and 

(6.24) 



Av 
2H n 



(R + AvS)p n (;v + )q n (;v-) -(R - AvS)p n (;v^)q n (;v + ), 



Av 
2H n 



= -(R + AvS)p n -i(; v+)q n -i(; V-) + (R - AvS)p n -i(;v-)q n -i(;v+), 



for n > and where the decoupling factor H n (u) has deg x H n — 0. Then the deformation coefficients of the 
deformed ID- semi- classical class are polynomials in x. 

Proof. We offer a proof in the spirit of that for Proposition 14.21 using Laguerre's method, and apply it to 
(|6~25)) first. From and ([5T2T1) we have 

=(R - AvS)f(; v+) - (R + AvS)f(; v.) - AvT, 
=(R _ AvS f n+P ^- 1 (; v+) - (R + AvS) qn+P< ^- 1 (; «_) - AvT, 

Pn Pn 

(R ~ AvS)p n (;v_)p^l 1 (;v + ) - (R + AvS)p n (;v + )p ( ^i 1 (;v^) - AvTp n (;v + )p n (;v_) 



(6.25) 



Pn(;V+)p n (;V-) 

(R - AvS)p„(; v-)q n (; v + ) - (R+ AvS)p n (; v + )q n (; u_) 



p n (;v+)p n (;v-) 

In this last expression the numerator of the first term is clearly a polynomial in x which implies the numerator 
of the second term is likewise. This latter numerator is an odd function under the exchange v + -h- v_ and we 
denote it by A«$„(i; u)/(2H n ), thus deriving (|6.23|) . The relation (|6.24l) is simply the case of (|6.23p under 
ti i — y 7i — 1. 

The remaining relations may now be derived from the following argument. From the workings of the 
previous paragraph we know that 



(6.26) 



AK> " (R- AvS)p n (;v-)p£± 1 (;v + ) - (R + AvS)p n (; v + )p { ^ «_) - AvTp n (; v+)p n (; «_). 



2H. n 



By multiplying the left-hand side of this relation by the Casoratian a n {v+)[p n -ip^}_ 1 — PnP^n—2\{'^ v +) = 1 we 
observe that the resulting equation separates into two terms with the factorisation 

(6.27) Pn (;v + ) 

= p£l 1 (;v + ) 



-a n ( W+ )^ P i 1 2 2 (; V+ ) + (i? + A^)pi 1 2 1 (;^) + A w r Pn (; U _) 



-a ra (t; + ) ^fc pn-i(;i> + ) + (R - AvS)p n (;V-) 



2H„ 



which implies that this expression contains the polynomial factors p n (;u+)p„_ 1 (;u+) and therefore can be 

written as p n (; w+)p^ 1 2 1 (; v+)tti, where m is a polynomial in x. Given that the two factors p n (; v + ), p| l 1 2 1 (; w+) 
are non-zero this then leads to two evaluations for tti , which by constructing a suitable combination of these 
and employing the Casoratian once more we find, upon some simplification, 

(6.28) 7Ti = a n (v + ) [-(R + AvS)p n -i(; v+)q n (; v-) + (R- AvS)p n (; u_)g n _i(; v + )} . 

Alternatively we could have multiplied the left-hand side of (|6.26j) by the Casoratian a n (v-)[p n -iPnli — 

PnP^n-2\{'i v -) = 1 an d deduced the factorisation 



3.29) 



-a n (v-) 



nJl) 



2H n 



Pn- 2 (; V-)-(R- AvS)p^l 1 (; v+) + AvTp n (; v+) 



p£i(;«-) 



-On(«-) 



Av$ r 

2H n 



■Pn-i(;v_) - (i? + Au5 , )p„(;i; + ) 



This fact then allows us to conclude that either of these expressions can be written as p n (; V-)p^_i{\ v-)ft2, 
where TT2 is another polynomial in x. Using an identical procedure to that employed above we can infer that 

(6.30) 7T 2 = a„(u_) [-(R + AvS)p n (;v + )q n -x(;V-) + (R- AvS)p n -i{j "-)<?«(; «+)] ■ 
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Relations (|6.20j) and (|6.21l) then follow from the definitions 



2R n - AvT n + 2H n - 



R + AvS 



2H n . nTT R-AvS 2H n 

-TTi, 2R n + AvT n + 2H n — — = — rAr^a, 



a n (v-j a„{v + ) a n (v + ) a n [V-) 

whereas (|rT2"2"j) follows from Av(T n + E n )/2H n = (R + AvS)/a n (v^) - {R - AvS)/a„{v + ). □ 

Analogous to the large x expansions of the spectral coefficients we have the following expansions for the 
deformation coefficients. 

Proposition 6.5. Let L = max(deg 2 ,i?, deg^S*). As x — > oo we have the leading orders of the terminating 
expansions of the deformation coefficients 



.32) — Rn = -(7„(«+) + 7»(0) 



H, 



R-AvS R + AvS 



7n-l(0 7n-l(0 



n-1 



£(&i(o-MO) 



i=0 



( i? + A^)-^±l T -(i?-A^)- 7 " (u - ) 



Aw 

^ 2^ = 



(i? + At-^)^^ - (i? - AwS) 



7n-l(0 

7n(0 



7n-l(0 



x" 1 + O^- 2 ) n>0, 



(i? + AvS) 



ln{v+) 

In(V-) 



7n(0 
71 n — 1 

5>(«_)-5>(« + ) 



i=Q 



i=0 



7n(0 



(R - AvS) 



7n(Q 
7n(0 



i=0 



i=0 



+ 0(a/'-' i ), n>0, 



Aw 

2H^ 



{R + AvS) 7ltzlM _ {R _ AvS fn-l(V-) 



(R + AvS) 



7«-l(Q 
7n-l(0 



7n-l(0 
n—1 n—2 

2 mo- E mo 



7n-l(0 



i=0 



i=0 



(i? - AvS) 



7n-l(0 
7n-l(0 



n-2 n-1 

$>(0-£&i(0 

,i=0 i=0 

+ 0(.x L ~ 3 ), n>0, 



Aw 

•35) — r„ = (7„(u-) - 7n(0) 



and 
(6.36) 



i? + AwS i? - AwS 



7n-i(v-) 7 n _i(u+) 



£(MO "MO) 



i=0 



(J? + AwS) 



7«(Q 

7„_i(u_) 



(i? - AwS) 



Aw, 

Hn 



(7„(u-) + 7„(0) 



i? + AvS R - AwS 



7n-l(«-) 7n-l(0 



5>i(v + )-&i(0) 



(i2 + AwS) 



7n(Q 
7n-l(0 



7n(0 



7n-l(0 



a;" 1 + 0(a; L - 2 ), n > 0, 



(R - AvS) 



7n(0 



7n-l(0 



a;" 1 +0(a; L - 2 ), n > 0. 



Proof. These formulae follow from the substitution of the expansions (|3.34[) and (|3.35p into the definitions 
(I6.20H6.24D . □ 

Remark 6.1. Unlike the case of the spectral coefficients (see Proposition I4.3[) the leading terms of the 
deformation coefficients are not determined by the weight data alone but depend non-trivially upon the 
three-term recurrence coefficients. 

Corollary 6.1. In the regular, deformed TB-semi-classical case the deformation coefficients R n , T„, $„, S n 
are polynomials in the spectral variable x with degrees independent of n, Aeg x R n = deg^Tn = deg,j.S n = 
max(deg a .iZ, deg x S) and deg^n = deg^^ = max(deg 2 .i?, deg^S*) — 1. 
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The expressions for the deformation coefficients (| 6. 2 411 6. 2 4")) can be inverted yielding a system of linear 
divided-difference equations in the deformation variable for the polynomials and associated functions. 

Proposition 6.6. The OPS corresponding to a deformed D- semi- classical weight satisfies the deformation 
divided-difference equation 

(6.37) O u Y n := B n M u Y n = i- (T* ®A U u Y n , n>0. 

-fin V w n ^n) 



Proof. The essence of this result involves the inversion of Proposition 16.41 which is carried out in a manner 
analogous to the proof of Proposition 14.41 □ 



Of the four coefficients r n , <£>„, S„ only two are independent because of the following relations. 
Proposition 6.7. The deformation coefficients satisfy the linear identity 
(6.38) <&„ = — —$„_i, »>1, 

0>n— 1 

and the trace identity 

'R + AvS R-AvS 



R- 


hAvS 




R - AvS 


a (v- 


-ho(v- 


) 


a (v+)'yo(v+) 


R- 


^AvS 




R - AvS 


a (v- 


-ho(v- 


) 


a (v + )-y (v+) 


R- 


\-AvS 




R-AvS 




-)7o(«- 


) 


ao(v+)lo(v+) 



(6.39) At- (r„ + H n ) = 2H n — — , n > 0. 

a n {v-) a n (v + ) 

Proof. The first relation follows by comparison of (I6.24[) with (I6.23|) and setting H n /H n _i = a n ja n -\. The 
second identity is easily seen from the expressions (|6 . 2 1 1) and (|6.22[) and has already been observed in the 
conclusions of the workings for the proof of Proposition ^. 41 □ 

The initial values of the deformation coefficients are 

(6.40) i? = -iffo[7o(^)+7o(«+)] 

(6-41) r = -L Ho { lo (v_)- lo (v + )} 

Av 

(6.42) S = -±-Hokyo(v-) + 7o(«+)] 

(6.43) $-i=0, $ = -2H 0l0 (v-)j (v + )T, 

(6.44) * = 0. 

Compatibility between the matrix recurrence relation and the deformation divided-difference equation 
implies the next result. 

Proposition 6.8. The recurrence matrix and the deformation matrix satisfy 

(6.45) B n+1 ■ M u K n - M u K n ■ B n = B u K n - \Av 2 B n+1 ■ D u K n ■ B n , n>0, 
or equivalently 

(6.46) K n (;v+) (l - ^AvB^ 1 (l + \AvB n ) = (l - ± At>5„ +1 ) ~' (l + \AvB n+x ) K n (;v-), n>0, 
Proof. The first form of the compatibility relation is deduced by comparing 

B u Y n+1 = B n+1 M u Y n+1 

= B n+1 M u (K n Y n ) 

= B n+1 (M u K n ■ M u Y n + ±Av 2 B u K n • DJ„) 
= B n+1 (M u K n + \Av 2 D u K n ■ B n ) M u Y n , 

and 

(6.47) B u Y n+1 = D u K n Y n = (D u K n ■ M u Y n + M u K n ■ B u Y n ) = (B u K n + M u K n ■ B n ) M u Y n . 

The second form is found by computing Y n+ i(x;v + ) in terms of Y n (x;V-) in the two possible ways corre- 
sponding to the orders of the operations jii->ti+1 and v— t— > v + . □ 
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Compatibility of the three term recurrence relation and the deformation divided-difference equation im- 
plies the following result which can be found by considering the representations of the deformation coefficients 
given in Proposition 16.41 

Proposition 6.9. The deformation coefficients i? n ,r n ,$„ satisfy recurrence relations in n, 
la ao\ a n+l(V-) f 0R , a -T ^ i an ( V ~h ou > a -p \ 

-(— 2-Rn+x + AvT n+ i) H — (2R n + AvT n ) 



TT V "T-i ' I'Tti ' tt 

tl n +l tl n 



, , ..Au, „ . . fR + AvS R-AvS\ 
= -[x- b n {v-)}— + 2a„(v_) ? — ? — - , n > 0, 



(6.49) a "+iK) (2i4+i + Awrn+i) + ^pl { _ 2Rn + AvTn) 

, , m A v^ „ , ,/R + AvS R-AvS\ 

= -[x - b n (v+) — $„ + 2a„(u+) , — r , — — , n > 0, 

V M u -J a n {v + ) ) 

Recurrences Ij6.48\ l and l{6.49{ ! are the analogs of combinations of l{4.43\ ! o,nd \4-44ty - 

Proof. The simplest proof of these relations is the one employing the definitions (|6.20I6.21I6.23[) of Proposi- 
tion 16.41 Taking the first relation we deduce 



o^v-) { _ 2Rn+i + AvTn+i) + ^l {2Rn + AvTn) + [x _ bn(i; _ ) ]^ $n 

tln+l tl n ti n 

= 2(R + AvS) [1 - p n (;v + )(a n+1 (v-)q n+1 (; w_) + b n (v- )?„(«_) + a n (v-)q n -i(v-) - xq n (v-))] 
a n (v-) 



- 2(R - AvS) 



2a n (v-) 



a n {v + ) 
R + AvS R-AvS 



qn(;v + )(a n+1 (v_)p n+1 (;v_) + &„(«_)p n (w_) + a„(w„)p n _ 1 (v_) -a;p„(w_)) 



a„(u_) a„(w+) 

where essential use of the three-term recurrence relations (13.513. 22[) has been made. The second relation can 
be found in an identical manner. □ 

We find bilinear identities for the deformation matrix analogous to those of Proposition 14.71 

Proposition 6.10. The deformation coefficients satisfy the bilinear or determinantal identity 

(6.50) Rl + \Av* [T n E n - = -H n R n \ R + A \ S + R ~ AvS ] , n > 0. 

L a n(v-) a n {v+) J 

Proof. A direct method of proof is possible substituting the expressions for the deformation coefficients 
Rn, r„, E n , $„, \f n in terms of products of polynomials and associated functions as given by (|6.20ti6.24j) into 
the left-hand side of the (|6.50j) . After expansion and considerable cancellation we recognise the the form for 
the right-hand side. □ 



From the above result we can deduce matrix identities analogous to Corollary 147 
Corollary 6.2. The matrix factor appearing in \6.4(fy has the determinant evaluation 

(6.51) det(l ± \AvB n ) = - 2Hn } R T . A f> , n > 0, 

R n a n {v±) 

and its inverse is 

la k<*\ rt j- i a„,b \-i a n {v±) ( Rn ± \Av E n T^Au $„ \ 

(6.52) (1 ± TrAvBn) = — ,„ t — ttt i a T r. 7 i a i n> 0. 

{ J V 2 ' 2H n {R^AvS)\ =f| Au *« R n ± ±Au T n J ' 

Proof. We note from the matrix formula for £>„ (|6.37p that 

det(l ± \AvB n ) = ^ {i? 2 + ±Aw 2 [r„S„ - ± \Av R n {T n + 3»)} 

2H n R =F Aw,? 

i?n O n (u±) 
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where we have used (|6.50[) and (|6.39[) in the last step. The inversion formula follows from this last result 
and the standard formula for an inverse. □ 



Consequent to the results of Corollary 16.21 we have the following expression for the matrix product ap- 
pearing in (|G.46I) . 

Corollary 6.3. The matrix product l[6.46\ ) has the evaluation 



(6.53) (l-±AvB n ) 1 (l + $AvB n ) 

( 



a n (v-) 



2H n {R + Av S) 



V 



2R n + 2H n R Al ' S + Av r n 
a n {v + ) 

Av V„, 



Av $ n 

2R n + 2H n R + AvS -Ac r„ 



n > 



\Pn-l\> v+) I 



a n (v-) 
R- AvS ( q n (;v+) 



R+AvS \q 



f.Z \ )®(Pn-l('>V-), -Pn(;v- 
l{i v +)J 



Proof. Using the inverse (|6.52p we form the matrix product and employ (|6.50p to simplify the diagonal 
elements, with the result given by (16.53[) . □ 



This result also motivates the following definitions 

(6.54) m+ := 2R n + 2H n R ~. Av f + Av T n , 

(6.55) «ft_ := 2R n + 2H n - + Av f — Av T n , 

a n (v-) 

(6.56) <£+ := -Av $„, := Av ¥ n , 
valid for n > 1, and for n = we have 

(6.57) D\+{n = 0) := -2 H \ lo( f + \ {R + AvS), 

ao{V-) 7o (y-) 

(6.58) «H_(n = 0) := -2 ~f° . ^°y~{ (R - AvS), 

ao{v+) 70 (v + ) 

(6.59) ¥+(n = 0) := 2A«tf o7o ( U+ ) 7o ( U _)T, <£_(n = 0) := 0, 
together with 



(6.60) S*.- ( ^ 
In these variables the deformation bi-linear relation (|6.50[) takes the form 

(6.61) dctS* = m+m- = , 4 f". N (i? 2 - At; 2 ,? 2 ). 

a n (u+)a„(v_) 

It is of advantage to write out recurrence-deformation compatibility equations in terms of 9^±,n,*P±,n, 
where again we append a subscript to indicate the dependence on the index n. We will henceforth denote 
variables evaluated at the advanced and retarded deformation variable by a n — a n (v + ),a n — a n (v-). 
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Corollary 6.4. Solving for B* +1 we deduce 

(6.62) 



1 



a„ + ia„ + i*P+,,i+i = -a n (x - 6„)9^ +! „ + a n (x - 6 n )9t-,n + (x - &n)(a: - &n)^P+,n + anting-, 



( 6 . 6 3) ^2_*p_ jn+1 = qj + „ ; 

(6.64) ^-W+.n+i = + ^^ +) „, 

(6.65) #^-,„+i = ^«+,„ - ^^+,„. 

The compatibility relation between the spectral and deformation divided-difference equations takes the 
following form. 

Proposition 6.11. The spectral matrix A n {x; u) and the deformation matrix B n (x; u) satisfy the Ui-Schlesinger 
equation 



(6.66) B u A n - B x B n + M u A n ■ M x B n - M x B n ■ M u A r , 



= \Av 2 U x B n ■ B u A n • (1 - ±Ay 2 Av 2 B x B n ■ B u A n ) 1 (M x B n + \Ay 2 D x B n ■ M u A n ) 

- iAy 2 M„^„ ■ B x B n • (1 - jgAy 2 Av 2 H u A n ■ B X B„) _1 (M u A n + ±Au 2 B„A„ ■ M x B n ) , 
or the equivalent form 

(6.67) (1 - \AyA n {;v+))~ 1 (l + \AyA n {- v+)) (l - |A«fl„(y_;)) -1 (l + \AvB n {y_- )) 

= (1 - \AvB n {y + ;))- X (l + \AvB n (y + ;)) (l - ±AyA n (; w-)) _1 (l + \AyA n {-v.)) . 

Proof. The first form follows from 3 x V) u Y n = 3JD) x Y n and the Leibniz formulae of (|2 . 10[) and (|2.1ip . from 
which we deduce 

[1 - ±Ay 2 Av 2 lB u A n ■ B x B n ] B X M U Y = [M u A n + \Av 2 D u A n ■ M x B n ) M U M X Y, 
[1 - ±Ay 2 Av 2 B x B n ■ B u A n ] B U M X Y = [M x B n + \Ay 2 D x B n ■ M u A n ] M X M U Y. 
Using equivalent forms of the spectral and deformation divided-difference equations (|4.35p and (|6.3Z[) we see 

Y(y + ; v+) = (1 - \AyA n {; V+))' 1 (l + ±AyA n (; v+)) Y(y_;v+) 

= (1 - \AyA n (-v + )Y 1 (1 + \AyA n (-v + )) (l - \AvB n {y_; (l + lAvB n {y^,)) Y(y_; V -), 

whereas 

Y(y + ;v+) = (1 - \AvB n {y + -)Y 1 (l + ±A W 5„(y + ;)) Y(y + ; V -) 

= (1 - \AvB n {y + - )Y 1 (1 + \AvB n {y + - )) (l - §AyA„(; _1 (l + \AyA n {- V -)) Y{y^v.). 

This gives us the form (|6.67l) . □ 

Remark 6.2. Let us denote the two fixed points of the x-lattice by xl and xr. By analogy with the linear 
lattices let us conjecture the existence of fundamental solutions to (|4.35[) about x — xl, xr which we denote 
by Yl , Yr respectively. Furthermore let us define the connection matrix 

(6.68) P(x;u) := (Y R (x; u))- 1 Y L {x- u). 

From (I4.35P it is clear that P is a D-constant function with respect to x, that is to say 

P(y+;u) = P(y-;u). 

In addition it is clear from (16.37)) that this type of deformation is also a connection preserving deformation 
in the sense that 

(6.69) P{x;v+) = P(x;v-). 

34 



The compatibility relation (|6.67p xB^iy+'i u )-^n( x l W-) = A*^{x] w+)i?*(y_; u) may be written component- 



wise in 


the new vBxicibles in the 


more prEictiCcil form 


c.lS 














(6.70) 


x[W + (x;V-)9\ + (y + ;u] 


)-T_(.T;«_) ( P + (y+ 




= 2U+(a:;t; + )9t + (j/- 


;u) 


-T+ 


(a;;w+)«P_ 


(y- 


\ u . 


(6.71) 


x[Z+(x;v-)m + (y+;u) 


+ W-(x;v_)¥ + (y+ 


;«)] 


= T+(a:;t; + )9t_(l'- 


;u)-\ 


-2B + 


.(x; 


-(y- 


;u 


(6.72) 




+ W+{x;v_)¥_(y + 


; w )] 


= T_(a;; U+ )SH + ( 2 /_ 




-2U_ 


(x;v+)¥^ 


(y- 


;u 


(6.73) 




)-1 + (x;v-)¥-(y+ 


;«)] 






-2L 


(x;v+)^ 


(y- 


;u 



where x 1S defined by (I6.17[) . The real content of (| 6. 7011 6. 73")) lie in the case n > 1 as the n = evaluations 
of (I6.72j) is trivially satisfied, whereas that of (|6.70|) and (|6.73l) are identically satisfied due to the definition 
(|6.17[) . and (|6.71|) is equivalent to the consistency identity (16. 16|) . 

Having assembled all of the ingredients of our theory we have to draw them together and to perform 
three tasks. The first is to parameterise the spectral coefficients in a minimal way consistent with the 
constraints of the spectral structures, namely (|4.46|) . The second task is to close this system, that is to 
say relate the deformation coefficients to the spectral coefficients, preferably with the parameterisation 
found in the first step. The third task is to derive the dynamical equations for this parameterisation with 
respect to the deformation variable. There are two avenues of approach to the problem of extracting useful 
information from the above compatibility relations (|6.70H6.73| . One way is to clear the denominators on 
both sides and resolve the resulting expressions in terms of © fc>0 x k + Ay © fe >o x k . This is useful for certain 
results, as we shall see in the application. The other method is to work with the evaluation of the spectral 
and deformation coefficients at certain key ordinates and construct a parameterisation based upon these 
variables. A fundamental and crucial role will be played by the zeros of the polynomials W 2 — Ay 2 V 2 , O ra 
and x as these ordinates, however these cannot be interpreted as singularities in the spectral plane. 

Remark 6.3. The compatibility relations (| 6. 70116. 75)) are satisfied identically at the common zeros of (W 2 — 
Ay 2 V 2 )(x;v±) in the sense that both sides of the relations are identically zero, that is to say at those fixed 
zeros independent of u. In addition one can show that the left-hand side of the bilinear relation (|6.61j) vanishes 
at the movable zeros in either of the two above cases. This implies that when (W 2 — Ay 2 V 2 )(x; v±) = 
then (R 2 - Av 2 S 2 )(y±;u) = 0. 

7. M — 3,i = 1 CASE AND DEFORMATION OF THE ASKEY-WlLSON OPS 

The previous sections treated the B-semi-classical orthogonal polynomial system with general divided- 
difference operators D, M and arbitrary degrees M, L for the spectral and deformation coefficients respec- 
tively. Here we apply the forgoing theory to the symmetrised form of the g-quadratic lattice in the elliptic 
sub-case, which singles out the Askey- Wilson weight and the non-trivial examples of deformations of this 
weight. We will consider the situation of the deformation variable also on this sub-case of the g-quadratic 
lattice. The deformation variable is u = (t + 1~ 1 )/2 with t = and the analogous relations for the nodes 
on the deformation lattice are 

(7.1) v++v_ = {q 1 / 2 + q- 1 / 2 )^ 

(7.2) v+-V- = W*-q- 1 '*)(t-t- 1 ), 

(7.3) Av 2 = {q 1 l 2 -q- 1 l 2 ) 2 {u 2 -l) 1 

(7.4) w+u _ =u 2 + l (g l/2_ g -l/2 )2 ^ 

7.1. Moments and Integrals. Our starting point is the simplest extension of the M = 2 case 

6 

(7.5) W ± AyV = JJ(1 - a,g- 1/2 z ±1 ), 

i=i 

which is a minimal, natural extension of the Askey- Wilson case given by (|5.7[) . The reader should appreciate 
that all of the conclusions we are going to draw will follow entirely from (|7.5p . so that any valid solution of 
the Pearson equation (14.21) for the weight with this data is as acceptable as any other. The most important 
consideration in selecting a solution is the support for the weight and existence of the moment data. As 

35 



we will see (i.e. in Propositions 17.21 and I7.7[) the only place where a specific choice for the weight enters is 
in determining the initial values of recurrences or divided-difference equations. The four fixed parameters 
ai, . . . , 04 appear in the same form as they do in the Askey- Wilson weight (M = 2) and so constitutes the 
base expression, and we seek to introduce an extra parameter and its associated deformation variable. We 
set 05 = at, ag = at~ x and designate a,t as the deformation parameter and variable respectively. 

Remark 7.1. One might entertain the possibility of a weight with the spectral data 

(7.6) W ± AyV = z^ 2 f[(l ^V^X 1 - q- l 'Hz ±l ){\ - g^ 1 / 2 ^ 1 ), 

3=1 

which is a M = 3 case on a q-quadratic spectral lattice, however it is q-linear in either of the deformation 
variables t, q a . This is in fact the spectral data for the weight in the integral representation of the very- well- 
poised %W*i basic hypergeometric function, see Eq. (1.13) of [5T]. This function is the classical seed solution 
for the Ej g-Painleve system as constructed by [52"] . 

Our first undertaking is to derive linear recurrences for the moments which will be required in the verifi- 
cation of key identities from Sections 2] and G3 to make contact with the n = seed solutions of the evolution 
equations derived later and for checking purposes. Through suitable choices of parameters k, b, a the moment 
recurrences can be recast as q- difference equations for an integral 73(01, . . . , ae) which is a generalisation of 
the Askey- Wilson integral (|5.24[) . 

Proposition 7.1 ([97 ). Let o~k denote the k-th elementary symmetric polynomial in a\, ... ,124 or a%, . . . , ag 

depending on the context. The integral I3 (a\, . . . ,ae) satisfies a three-term recurrence in one variable, which 
we take without loss of generality with respect to ax, 

(7.7) Q = JJ(aia i -l)/ 3 (oi J ...) 

6 5 

+ [1 + q~ 1 - aiC^aj - qa!) +Y[a j (a 1 ^2a~ k 1 - q^ 1 - (q + l)a?)]7 3 (gai, . . .) 
3V1 1 kjti 
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1 

In addition the integral 7 3 (ai, . . . , <i6) satisfies a three-term recurrence in two variables, taken to be with 
respect to 05, oq, which constitutes a pure recurrence in the deformation variable u 

4 

(7.8) = (a 5 - qa 6 ) J|(l - a J a 6 )7 3 (. . . ,q 2 a 5 ,a 6 ) 
3=1 

- (05 - a 6 ) [(1 + g)(l + qa 5 a 6 a 2 + q 2 a\a\o±) - (qa 5 - a e )(qa 6 - a 5 )(q + 04) 
-q(a 5 + a 6 )(<ri + qa 5 a 6 a 3 )} 7 3 (. . . , qa 5 , qa 6 ) 

4 

+ (qa 5 - a 6 ) JJ(1 - aja 5 )I 3 (. . .,a 5 ,q 2 a 6 ). 

3 = 1 

Proof. The first recurrence can be read off from (|4.11[) specialised to k = 0, after evaluating the expansion 
coefficients 

x i \ a6 ~ 1 
03,2 (a) = . , 
q°' z a z 

^3,i( a ) = s/o A Qafa - ci) - (1 + q)( l + ?a 2 )( cr 6 - 1)], 
q°' z a z 

*3,o(o) = "wirK 1 + fl2 + a4 )( CT 6 - 1) - (a + a 3 )(^5 - cri) + a 2 (o- 4 - cr 2 )]. 
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The second recurrence follows from combining the first with the general identity 

a k I(qa v ...,a k )- ajl(aj, . . .,qa k ) = (a k - aj)(l - a i a fc )/(a i , . . . ,a k ), 

which applies to any integral with products of </>oo(^; «j)</'oo(^; o-k) m the denominator of the integrand and 
any distinct pair cij 7^ a k . □ 

Remark 7.2. A particular solution of (|7.8p which serves as a concrete example of a moment sequence is one 
taken from [57] 

j, 2 -1,-1 1/2 , 1/2 -1.-1 \ {a 1 1 a i at,a 2 1 (j i at,a 3 1 a i at,a i 1 aiat\q) 00 

(7.9) m 0fi (t) = t 1 (qat,a t ,q 1 at,q 1 a t ;q) 00 — — 3 „ 5-5— 

[axca, a2at, a 3 at, a^at, t , 174 art , q)oo 

x sW7(q~ 1 a4a 2 t 2 ; q~ 1 <i4a 2 , a\at, a^cd, a 3 at, a^at; qa~ 2 ) + (t i-> t^ 1 ), 

although this is not a general solution of the q-difference equation, given subsequently as (|7.78p . 

Henceforth we denote dj as the j-th elementary symmetric function of the parameters g _1 / 2 ai, . . . , q~ 1 ^ 2 aQ. 
If <r k denotes the fc-th elementary symmetric function of a\, a%, 03, 04 then 



(7.10) 


cti 


= g" 1/2 (cri + 2cm), 


(7.11) 




= q~ 1 ((T2 + a 2 + 2aa\u), 


(7.12) 


5-3 


= q~ 3/2 {a 3 + a 2 oi + 2aa 2 u), 


(7.13) 


(T 4 


= q~ 2 ((T4 + a 2 a 2 + 2ckt 3 u), 


(7.14) 


CT5 


= q~ 5/2 (a 2 a 3 + 2aaiu), 


(7.15) 


<7(i 


= q a (74. 



We will find it advantageous to define two "analogs" of the integers or half-integers s by [s] := q s a§ — q~ s 
and {s} := q s <Js + q~ s . For convenience we employ the notations w± for the evaluations w± = W(±l) = 
±1 — cti ± (T2 — ct 3 ± a 4 — (T5 ± We also require the coordinates Xj and defined in terms of the 
parameters by xj = \(a.j + aj 1 ) and ij — |(g _1 ' 2 aj + q^-^aj 1 ). Finally let us define the polynomial 
w(z) := Ilj=i(l — q~ 1 / 2 ajz) which is not to be confused with the weight w(x; u). 

Definition 7.1. Note that we have six free parameters at our disposal, 01, • • • , 014, ot € C, n € Z>o and one 
variable t subject henceforth to the following generic conditions - 

(i) q * 1, 

(ii) \q~ 1/2 aj\ ^ 1 so that Xj $ (-1, 1) for j = 1,.. . ,4, 

(iii) \q- 1 / 2 a\ ^ 1 so that g-V 2 a + g 1 / 2 ^ 1 ^ (-1, 1), 

(iv) a^0,±gV2 

(v) t ^ ±1, 

(vi) [n],[n+i]^0. 

Situations where one or more of the above conditions are violated have to be treated separately, which we 
refrain from doing here. 



As a consequence of 1|7.5[) we have 

(7.16) W[x) = 4(1 + & 6 )x 3 - 2(cti + vr )x 2 + (a 2 + <t 4 - 3 - 3<7 6 )x + a x - a 3 + a 5 , 

(7.17) V(x) = ^1/2 h 4 (l ~ ^e)^ 2 + 2(*i - <7 5 )x + 1 - ct 2 + ct 4 - ct 6 ] , 
The theory detailed in Section 0] allow us to evaluate the third spectral data polynomial. 
Proposition 7.2 ([97]). The polynomial U(x) is given by 

(7.18) U(x) = - g-1/2) [ m °'° ("^ + 91/2,55 - ^ 1/2 ^) - W( TO o.+ + m o,-) 
where mo,± = /Diw(i;m)2 ± . 
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Proof. We start with 

(7.19) AyU(x) = (W — AyV)f(y+) - (W + AyV)f(y-), 

and employ the proof in [97] . Having achieved the task of expanding the right-hand side of (|7.19p in canonical 
basis polynomials, all we require are explicit expressions for the coefficients given below 

(7.20) K 3i2 (a) ' 

(7.21) K 3i3 (a) 
The result is ((77151) . □ 



2g 9 / 2 a 3 
l + Q-6 
2q 9 / 2 a 3 ' 



[(i + « + r)(i + f«')(l + ffe) ~ ro(ffi + o*)], 



7.2. Spectral Structure. Our next task is to construct a minimal parameterisation of the spectral matrix 
A n (x; u) and is the subject of the following proposition. We also recall the large x expansions for the spectral 
data and coefficients of a more explicit nature which will be of importance 



(7.22) 
(7.23) 

(7.24) 
(7.25) 



W ± AyV ■ 



w±(q 1 / 2 -q- 1 / 2 )v 



(1 + q n )(l + q- n )w + (q 1 ' 2 - q- 1/2 )(q n - q- n )v 



x M + 



n„ + v 



qn+l/2 _ -n-l/2 

~~q 1 ^~~Trr- 

q n _ q -n 



U ; + (g»+l/2 + g -n-l/2 )u 



„M-2 



l/2_ q -l/2 W + (g ) 



„M-2 



gV" - g 

for some constants w, v independent of n,x. 

Proposition 7.3. Let one of our free variables, X n (u), be the zero of Q n (x;u) with respect to x. Also define 
the other independent variables for n > by 

(7.26) v n (u) = 2W n {X n ;u) - W(X n ;u), fi n (u) = Cl n (X n ;u) +V(X n ;u), 
which are related by the spectral conic equation 

(7.27) v 2 n - W 2 {X n ) = A(X 2 n - 1) [/4 - V 2 (X n )} , 

again valid for n > 0. Let us assume that X n {u) =^ ±1 for all n, u. The spectral coefficients for the M = 3 
deformed Askey-Wilson OP system are parameterised by X n , and either v n or fj, n , through the expressions 



(7.28) 2W n (x;u) -W(x;u) 



x 2 -1 



V n + (x- X n ) 



4{n}(x 2 - 1) + ± w+ -^— + l w _^—— 



(7.29) Q„(x; u) + V(x; u) = /i„ + 4 /2 W g _ 1/2 a?(a - A„) + ^ ^ 6 ff _i/ 2 — 



+ 



{n} 



gV2 - g-l/2 



g [n 



CTlO- 6 + CT 5 



A 2 -1 



4ct 6 

1 1 
1 — W- 

1 - A„ 2 1 + A„ 



(x - X n ) 
1 



and 
(7.30) 
for n > 0. 



8„(af; u) = 8 /2 _ g 2 _i /2 - An), 



Proof. Firstly we recall the spectral data polynomials given by (|7.51 17.1617.171) . The product W 2 — Ay 2 V 2 
plays a significant role and therefore we define another set of elementary symmetric polynomials by 



(7.31) W 2 - Ay 2 V 2 = Coo - Xj) = [x 6 - ax 5 + e 2 x 4 - e 3 x 3 + e 4 x 2 - e b x + e 6 ] , 



where Coo — 64ct6. We parameterise the spectral coefficients in the following way 

(7.32) 2W n - W = w 3 x 3 + w 2 x 2 + w x x + w , 

(7.33) 6„ = zu + (x - An), 0„_i = m-{x - A„_i), 

(7.34) Q n + V = v 2 x 2 +vix + v a , 

where the leading order coefficients in each are trivial and given by using (|7.16[) and (|7.17p in the expansions 
(ITm g2H) and $m 

( 7 - 35 ) w 3 = Hn}, OT+ = 8 ^"*^ , «2=4 gl/2 j"^_ 1/2 . 

From the fundamental bi-linear relation (|4.46l) we get a system of quadratic polynomial equalities 

(7.36) wl-Avl = C 00 , 

(7.37) 2w 2 w 3 — 2AwiW2 = — Cc^ei, 

(7.38) 2wiw 3 + k; 2 — A(v x — «2 + 2v v 2 ) + a^An7 + n7_ = Cooe2, 

(7.39) 2wiw 2 + 2wqw 3 - 2A(v vi - V1V2) - a 2 n Azu + zu-{X n + A n _i) = -Cooe3, 

(7.40) wf + 2w w 2 - A(t;g — v 2 — 2vqv 2 ) + a 2 l Aw + W-(X n X n -i — 1) = Cooe4, 

(7.41) 2w W! + 2AwgWi + a 2 n Auj + Ta_(X n + A„_i) = -C^es, 

(7.42) w 2 + Av 2 - alAzu+zu-XnXn-i = Cooe 6 , 

where A is, again, defined as (q 1/2 - q^ 1 / 2 ) 2 . Now (|7.42I7.40I7.38|7.36|) imply (w + w 2 ) 2 + (wi + w 3 ) 2 = 
Coo(l + e 2 + e 4 + e 6 ) while (|7.41|7.39|7.37| l imply 2(w + w 2 )( Wl + w 3 ) = -Coo(ei + e 3 + e 5 ). Forming the 
sum and difference of these two later relations we conclude that 

,„. , eWYlJ + eWC-l) eW(l) - eW(-l) 

(7.43) Wl +w 3 = i w n +w 2 = i '-, 

where e,e' = ±1 and are yet to be determined. Using (|7.43|) along with the above definition we can solve 
for wo,w 2 , which only leaves the signs e, e' unresolved. These can be fixed by requiring that the n = 
evaluation of 2W n — W precisely reproduces W . This is identically true for all free parameters a\, . . . , ae 
provided e = +1, e' = — 1, and yields (I7.46I7.44|) . The above coefficients are given by 

(7-44) w 2 = -p— - w 3 X n - \ W{1)—^— + 1 , 

A n — 1 A n — 1 A„ -f- 1 

(7.45) Wl = - W3+ iw(l)-%W(-l), 

(7.46) w = + w 3 X n + 1^(1)^- + iW(-l)- A " 



At-l — 2" "X n -1 2 " v y A„ + l' 

We observe that this condition also gives vq = W(Xo), as it must. Through knowledge of w 2 and utilising 
(|7.37|) we can determine Vi, which is given in (|7.47l) . 

(7-47) *»vi = ^- w 3 X n + 1^1- \W{1)-^— + 1 W (-1)-^—, 

w 3 Xi - 1 w 3 z X n - 1 A„ + 1 

However to find the remaining coefficient vo we require and this result is given by (|7.48p . 

m AQ \ Av 2 Av 2 X n Coo , 2 1 Cooei i wf1 i A„ , A„ 

(7.48 W = ~ T2 7"n + K ~ 2 A„ + 2 W{1)- - ±W(-1) ■ 

w 3 w 3 Xi - 1 w 3 w 3 " A„ - 1 A„ + 1 

When we examine Q n + V at n = (recall that fio = 0), we nn d in addition to previously found relations 
the equality [jlq = V(Xa), again confirming our definition. Therefore we have succeeded in relating the 
sub-leading coefficients explicitly in terms of two independent variables. □ 

We now address the question of representations for the three-term recurrence coefficients. We shall find 
that another set of variables, although equivalent to v n and /i„, will lead to the simplest forms for the 
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relations we seek. 

(7-49) w 2 , n := - 4{n}A„ - ^(1)^ + 

(7.50) v , n ■■= (q 1/2 - q~ 1/2 )Hn - 4[n]A£ - A„ ({n}w 2 ,„ + 4(<7i<7 6 + a 5 )) . 

M 

To achieve this outcome we require the following recurrence relations. 

Proposition 7.4. T/ie dynamical variables w 2tTl {u) and i>o,n(u) satisfy the following system of first order 
coupled recurrence relations in n 

(7.51) w 2 , n +i - w 2 ,n = ~^{q 1/2 - q- l/2 )[n + |]A„, 

IVI + — 1 / \ 

(7.52) v , n +i + v ,n = - T — 2\ 2{n + \}\ n w 2 , n + 8[n + ±][n]A 2 + 4(g 1/2 + <T 1/2 )(aiCT 6 + ct 5 )A„ , 

[n + 1J [nj V - / 

and vaZid /or n > and subject to the initial conditions U(Xo;u) — 0, W2,o = — 2(<Ji + (75) and «o,o = 

1 — (72 + (T4 — (T6 . 

Proof. For the first recurrence we employ our expressions (|7.28[) and (I7.30|) in (|4.43f) and equate the coef- 
ficients. The coefficients of x 3 and x are identically satisfied whereas that of either x 2 or a; yields (|7.51|) . 
The second recurrence follows from the examination of the x terms in (|4.44j) and the employment of the 
first recurrence. We note the terms in x 2 cancel identically, as they must. □ 

Explicit evaluations of the three term recurrence coefficients can be given in terms of this alternative set 
of variables. 

Proposition 7.5. The three-term recurrence coefficients are found to be given by 

(7.53) 16[n + \][n] 2 [n - ±]a 2 = o % w\ n + 2(oict 6 + v 5 ){n}w 2 , n + 2[n] 3 w ,„ 

+ 2[n} 2 (4ct 6 + 2ct 4 + 2ct 2 ct 6 + 2a x a 5 - {n}(l + a 2 + & i + a 6 ) + 2[n} 2 ) + 4(ct 1 ct 6 + a 5 ) 2 , 
for n > assuming a 2 , = 0, and 

(7.54) [n + l)[n)b n = -\{n + \}w 2 ^ n -[n+ \][n]\ n - \{q 1/2 + q- 1 ' 2 )^^ + cr 5 ) 
again valid for n > 0. 

Proof. Using (I7.38P we can solve for a 2 . Observe that the right-hand side of this equation is independent of 
n whereas individual terms on the left-hand side are. If we assume a 2 , = then the right-hand side is equal 
to that of the left-hand side evaluated at n = 0, which we know because we can express it simply in terms 
of the parameters. This means a 2 is expressible as a sum of differences and after considerable simplification 
we arrive at the expression (|7.53[) . To find (|7.54p we start with (|4.44[) and use our previous results for (|7.30p 
and (|7.29[) . Examining the terms in x, we find an expression for b n in terms of X n ,y n and A n +i, fji+i) or 
equivalently in terms of \ n ,W2,n and A n +i, W2, n +i- Using (17.51[) we can eliminate the w 2 , n+ i term and after 
some factorisation we arrive at (|7.54p . □ 

We will not pursue the theory for the n 1— > n + 1 recurrences any further here but refer the reader to |96j . 
Consequently we will drop the n subscript from most variables until Subsection 17.41 

The forgoing parameterisation of our system given in Proposition [7j3] is not useful in the investigations of 
the u or t-evolution of our system and we require an alternative construction. In conformance with this we 
define the auxiliary variables l(t) and %±(t) by 

(7.55) A^Z + r 1 ), 

(7.56) 3 ± ^^liq^-q-^ll-r 1 }^ 

where, in the first case, the inversion is given by the branch whereby I — > 00 when A — > 00 and an identical 
choice is made for the second case. 
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Proposition 7.6. The spectral matrix elements have an alternative parameterisation 
(7.57) = 2a n [n + \]{z ~ z- 1 )^ ~ 1){1 - l^z' 1 ), 

and 

- z- x ){z - l- l ){q n u 6 - q~ n lz- 1 ) (z - z" 1 )^ - l)(q n a 6 - q^l^z' 1 ) 



(7.58) [n]W+(z;t) = j+- 



(I -I- 1 ) 2 °~ (I -I' 1 ) 2 

+ [n](z - z-^iz - 0(1 - r 1 ;*- 1 )^"^ - g _n z _1 ) 

+ [a 6 CTi + 5-5 - 2ct 6 (z + r x )](z - - i)(i - r 1 * -1 ) 



i (z + l)(z - 0(^~ 1 - Q(g"g 6 - q- n z- r ) i (z - l)(z - /)(z~ 1 - Q( g "g 6 + ^-"z- 1 ) 



(7.59) [ra]2D_ (z;t) =3- ^ _ pj^ +3+ » _ ^_ 1)2 

- [n](z - z- x )(z - 0(1 - r 1 z' 1 )(q n a 6 z- 1 - q~ n z) 

- [er 6 oPi + ct 5 - 2ct 6 (z + r x )](z - z-!)(z - 0(i - r 1 * -1 ) 

(z + l)(z-Q(z- 1 -0(g"q 6 z- 1 - g -") x (z - l)(z - Otz- 1 - 0(g"a 6 z- 1 + q- n ) 

Proof. These formulae can be viewed as Laurent interpolating polynomials satisfying the following evalua- 
tions at the given nodes 



(7.60) 


x±(z -1 ;*) = -%k(z;t), 


(7.61) 


1+(l ±1 ;t) = 0, 


(7.62) 


X+(±l;t)=0, 


(7.63) 


X+(z;f) - 2a„[n+±]z 2 , 

z— >oo * 


(7.64) 


%+(z;t) 2a n [n+\]z- 2 

z— >0 


and 




(7.65) 


W ± (z- 1 :t)=W T (z;t), 


(7.66) 


2H±(/;t)=3±, 


(7.67) 




(7.68) 


W+{z;t) ~ q n v 6 z 3 , 

z—^oo 


(7.69) 


W+(z;t) ~ n q~ n z- 3 . 

Z— ¥0 



The above formulae are also a consequence of the parameterisation (|7. 28117. 29|) . □ 

Remark 7.3. We note that the eigenvalues of A*(z;t) as z — ► 0,oo are q n a 6l q~ n whilst those at the fixed 
points of the lattice, of the matrix ^4*(±l;t), are w±. Under the mapping interchanging the interior and 
exterior of the unit circle in the spectral variable z H> 1/z we observe that 1± M> — 1± and 2U± n- 2H T . As 
a consequence A* n (z; )A* (z _1 ; ) = (VK 2 — Ay 2 U 2 )Id, and therefore the mapping corresponds to a reversal 
of direction on the spectral lattice. With respect to the mapping I H > I' 1 we note that this must be taken 
with 3 + and conclude that T + ,2tJ± are invariant under this type of transformation. In addition A* is 

symmetrical under 1 1— > £ . 

7.3. Deformation Structure. The deformation data polynomials R(x;u), S(x;u) are computed as 
(7.70) R± AvS= (l-q- 1 / 2 at ±1 z){l-q-^ 2 at ±1 z- 1 ) = fc&q-^at^ 1 ), 



which implies 

(7.71) R=l- c?q- x - 2aq- 1/2 xu + 2a 2 q- 1 u 2 , 

(7.72) S = 2^—{aq- 1 ' 2 u-x) 1 

q - 1 

indicating that L = 1. We observe that the polynomial 

(7.73) R 2 - Av 2 S 2 = 4g- 1 a 2 (^ - x 5 ){x - x e ), 
divides 

4 

(7.74) W 2 - Ay 2 V 2 = G4q- 3 a 4 a 2 ]J(x- Xj) ■ (x - x 5 )(x - x 6 ), 

3=1 

in conformity with Remark 16.31 The role of the common zeros X5 , xq will be crucial in the ensuing investi- 
gations. We also note that these points can be represented on a u-lattice thus 25 = E~x$ and x§ = E+x§. 
We can compute the remaining deformation data polynomial T(x; u), which will be of degree zero. 

Proposition 7.7 ( 97 ). The deformation data polynomial T(x;u) has the evaluation 

(7.75) T(x; u) = J?L tm °^-) Z g^oCgjO , 

q — 1 t — t 1 

Proof. We start from the formula 

(7.76) AvT(x; u) = (R- AvS)f(x; v+) - {R + AvS)f(x; 

and employ the expansion (I3.46|) for some parameter a. Setting a to one of 01,02,03,04 and utilising the 
identity 

4>i{x]b) _ (l-o6g n )(l-6o- 1 g-™) ba^q- 71 
<f> n+1 (x;a) (j> n+ i(x;a) (j) n (x;a)' 

we find, by equating coefficients, the result (|7.75l) . Note that m doesn't depend on a. If one chooses a 
from 05, 06 then one requires the additional identity 

1 _ (1- q n )(a 2 - q- n ) q~ n 
<j) n (x;qa) <j) n+ i(x;a) <j> n (x;a)' 

in order to merge the two series in (I7.76p . Alternatively one can compute T from the initial value (|6.43[) 
using the asymptotic expression for $0 a s given by (|6.33[) . □ 

As a check we can verify that the spectral data polynomials W, V, U and the deformation data polynomials 
R, S, T satisfy all of the consistency relations formulated in Section |6] 

Corollary 7.1. The spectral data polynomials W, V, U and the deformation data polynomials R, S,T satisfy 
the consistency relations \6.4\l and H6.16\) . 



Proof. The first task involving W, V and R, S is elementary, whereas the second requires a moment relation 

(7.77) [l](mo,+ + mo,_)(t) = [q~ 2 a 2 o- 3 - q~ l ai + g~ 3 (g - a 2 )^ 2 ^ 1 *" 1 + a 4 at)] m , (*) 

_ wjq^^at) q-^Hmo^g-H) - q^H^m^jt) 
at q-VH-qVH-i 

and the moment g-difference equation in t 

(7.78) (q-H - gt" 1 )w(ar 1 )TO 0i0 (g 1/2 t) 

- q- 1/2 {q- 1/2 t - q 1 ^- 1 ) [(1 + q)(l + q- 1 a 2 a 2 + q- 2 a 4 a 4 ) + q- 1 a 2 (q + <j 4 )(t - t" 1 )^ 1 * - qr 1 ) 

-a(q-^H + g 1 / 2 t- 1 )(a 1 + g^aVg)] m ,o(<r 1/2 t) 

+ (i - r 1 )u)(<? -1 at)Trao,o(Q' -3/2 *) = 0, 
which is just (|7.8|l with a 5 h-> q~ 3 / 2 at, a 6 1— > q~ x l 2 at~ x . □ 
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We know for L = 1 that the deformation coefficients can be parameterised thus 

(7.79) 9\± = r 1± x + r ±, 

(7.80) ¥±=p±- 

Furthermore we know from the large x expansions (|6. 32116. 36]) what the leading terms are, however in contrast 
to those of the spectral coefficients these are related to the three-term recurrence coefficients a„,7„ in a 
nontrivial way. 

Corollary 7.2. The leading order terms of the elements of the matrix J5* are given by 

tln(v+) 
7n-l(«-)' 
ln(v-) 

fry„_i(«+)' 



(7-81) 




= Aaq- 1 / 2 ^ 


(7.82) 


n- 


= iaq-^Hn 


(7.83) 


p+ 


= Aaq~ l ' 2 H n 


(7.84) 


P- 


-Aaq- 1 ' 2 ^ 1 



7„(u_) tj n (v+) 

tj n -i(v + ) 7„_i(«_) 



On-1 



7n-l(«-) *7n-l(«+) 



/or n > 1 whilst for n = one can wse 7_i = ao7o- 



Up to this point we haven't exercised a choice regarding the decoupling factor H n but one can take 
H n — \o^n henceforth. 

Remark 7.4. We note that the eigenvalues of -B* (z; t) as z — > 0, oo are ri±. By construction B* n is symmetrical 
with respect to z i— >• z _1 . Under the mapping interchanging the interior and exterior of the unit circle in the 
deformation variable t i-> 1/t we observe (see ()6.54|) . (|6.55|) and (|6.56|l ) that *p± i-> — *p± and i— > JH^. 
As a consequence 

(7.85) B^jt^Ot" 1 ) = 4 f". . (i? 2 - A^^)Id, 

and therefore this mapping corresponds to a reversal of direction on the deformation lattice. 

As part of an efficient methodology we will require formulae which relate the lower order terms, with 
respect to x, of the deformation coefficients to the leading order term. 

Proposition 7.8. The trailing terms of the deformation coefficients are related to each other and to the 
spectral parameterisation by 



(7.86) 


P+ 


= ri+a, 


(«_) - ri 


-a n (v+), 




(7.87) 


P- 


= n+a, 


(v+) - n 


-a n (v-), 




(7.88) 


ro+ 


n+ 

Av 2 


-w 2 (v + ) 


+ w 2 (v-) 


- 2VA(q- 1 / 2 a + q^a-^d^H - q^r 1 ) 


(7.89) 


ro- 


ri- 
Av 2 


w 2 (v + ) - 


102 (W-) - 


2VA(q- 1 / 2 a + q^a- 1 )^^- 1 - q- n t)] , 



Proof. The formulae given above are ones which involve both shifts up and down on the u-lattice of the 
spectral variables, and these are easily derived by resolving the compatibility relation (|6.67[) as a Laurent 
polynomial in z. Such a polynomial is constructed in our application by substituting the elements (|7.571 
I7.59|) and (|7.79I7.80[) into (|6.67[) and collecting terms. For the 1, 2 component we find that the leading order 
non-zero contributions occur at z ±r and the coefficients of both these terms will vanish if 86[) holds. In the 
case of the 2, 1 component the z ±r terms are the leading non-trivial ones and these both vanish when (I7.87[) 
holds. The two latter relations (|7.88I7.89[) follow from the requirement that the coefficients of z 7 terms of 
the 1, 1 and 2, 2 elements vanish respectively. Finally we observe that an independent way of verifying (17.86)) 
and (|7.87p is through a trivial combination of the formulae (|7.8H7.84")l □ 

A vital part of our analysis consists of resolving the compatibility relations (| 6. 70116. 73]) with respect to the 
zeros and poles of x(z; t). 
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Corollary 7.3. The residues of the compatibility relation given in Provosition 1 6. 1 1\ consist of the following 
equations: 

(i) At the advanced co-ordinate z = zs = at 
(7-90) n-E u x 5 + r - = -— T797TP+' 



(7.91) r 1+ E u x 5 + r 4 



(ii) at the advanced co-ordinate z = E 2+ zq — q~ 1 at~ 1 

{im) ^ X6 + = Z + (El+z^H) P+ > 

(7 - 93) ri+K X6 + r ° + = I_(£g+z 6;g i/2 f ) P -' 

(iii) at the retarded co-ordinate z — E 2 ~ z$ = q~ 1 at 

1+(Etz 5 ;q-VH) 



(7-94) r 1+ E~x 5 + r 0+ = 5 '"_ 1/9 ,> +, 



(iv) and at the retarded co-ordinate z — zq — at -1 
7.96 ^l+^u ^6 + 7-0+ = -=r- ( _ , . p+, 

(7-97) n-E+ x 6 + r - = ——. zrn^rrP-- 

j4Z/ the deformation parameters are evaluated at t,u, i.e. p±(t),r\±(t),rQ±(t). 
Proof. The multiplicative factor (|6.17[> is given, in this case, by 

(7 - 98) (l-q-iatz-i)(l-at-i Z y 

At the zeros of \ the right-hand sides of (16.7016. 7116. 7216.73)) must vanish whereas at the poles the left- 
hand factors must vanish. Thus for the zero z — z$ = at applied to (j6.70|) we deduce (17.901) . For the 
other zero at z = l/E 2+ ze = qa~ 1 t applied to (|6.70[) we make additional use of the symmetries of A* n and 
B* with respect to inversion of z to derive (|7.92p . At the pole z = E 2 ~z 5 = q~ 1 at the residue of (|6.70[) 
yields (I7.94j) whereas at the other pole z = 1/zq — a~~ 1 t we have to apply the additional symmetries to the 
residue to arrive at (|T.96|> . However we have further relations, which are equivalent to the four relations 
derived above, by specialising the spectral variable in these ways due to the fact that (W 2 — Ay 2 V 2 )(z — 
at, q-^-ar^lH) = 0, (W 2 - Ay 2 V 2 )(z = q~ l a% a*" 1 ; <T 1/2 *) = 0, (R 2 - Av 2 S 2 )(z = q-V 2 at; t) = and 
(R 2 - Av 2 S 2 )(z = q- 1 ^ 2 at- 1 ;t) = 0, which implies 

(7.99) (2B+2XJ- + T+T-)(z = at, q- y aT y \ q 1/2 t) = 0, 

(7.100) (2B+2H- +T+T_)(z = g _1 at,at _1 ;g- 1/2 t) = 0, 

(7.101) (<K+<K_ + qj+ ( P_)(z = q- 1/2 at; t) = 0, 

(7.102) (m+V\- + <p+<p_)(z = (jf-^at -1 ;*) = 0. 

This means that each of the identities (|7.90[ 17.921 17.941 17.961) can take four forms and the remaining relations 
(|7.911 17.931 17.951 17.97[) are examples of just one of those forms. This equivalence also ensures that if the 
residue condition is satisfied by (|6.70[) then it is automatically satisfied by (|6. 7116. 7216. 73[) as well. □ 
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Corollary 7.4. The following product relation holds amongst the spectral coefficients evaluated at the zeros 
of the spectral determinant 



(7.103) 



W + (z&qVH)VB-(E%rz5-,q-VH) _ 2P + (z 6 ; g" 1 / 2 *) 2P^ (El+z 6 ; q^H) 
Z+faqVH) % + (El~z 5 ;q-i/H) 1+(z e ;q-VH) 1+(E 2 U + z 6 ;qVH) ' 



Proof. By comparing the right-hand sides of the pairs, (jTTMl 17351) . (I7.92l7.97p . (EH EH]) and (f7TMl 17^31 . 

we see that 



(7.104) 
(7.105) 
(7.106) 
(7.107) 



2P+(^;g 1/2 t) 
T+(2 6 ;« 1 / a t) 

2P-(£?g+z 6 ;g 1/2 f) 
1 + (E* + z 6 ;qi/ 2 t)' 
W.{El-z b -q^/H) 
1+(El-z 5 ;q-VH) 
W+jze-g^/H) 



P+ = 



%-{E 2 u -z b ;q-^H) 
W-jz^q-VH) 
<Z^z 6 ;q-VH) P -' 

2P,(z 5 ;g 1/2 

W + (E*+z 6 ;qVH) 
H-iEl+z^qVH) 1 



V- 



By forming the cross products of the first and second, or equivalently the third and fourth we can establish 
()7.103[) . assuming p+p- ^ 0. Taking cross products of the first and third, or second and fourth leads to the 
identities that are trivial consequences of the determinantal relation. □ 



Definition 7.2. In conformity with the remark made immediately preceeding Corollary 16.41 let us define 
the variables with deformation arguments evaluated at the advanced and retarded co-ordinates through the 
notation l^l^q 1 / 2 ^, I = /(g^ 1 / 2 ^ etc. 



Lemma 7.1. One solution for the component p + is given by 



(7.108) q^a' 1 ^- 1 -t)[n + 



li ^l- Or, 
11 P+ 



qtl 



(I -I- 1 ) 2 



3+ 



3- 



(at - l)(qt - al) (-a + qtl)(-l + cdl) 
1 (r n +9 w+1 *Ve) + M 1 )- ' 



at 



A 



A + l 



A second solution for the component p + is 



(7.109) q L/z a- L (t- 



t)[r 



qtl 



P+ 



(I -I- 1 ) 2 



(g- atl)(-a + tl) (t - al)(-at + ql) 



- f (q-' l t 2 + g" + V 6 ) + \w(l)J— + H-1) T 
at A - 1 



A + l 



Proof. The first relation (|7. 108[) follows from the subtraction of (|7.92|) from (|7.90|) and simplifying. The 
second relation (|7.109[) follows from similar reasoning applied to (|7.94l) and (|7.96p . □ 

Remark 7.5. A profound difference between the expressions (|7. 10817. 1091) and the corresponding formulae 
in the case of the g-linear spectral lattice is the presence of both j± (or j±) whereas in the latter case only 
one these variables in present, which significantly complicates the ensuing analysis. 

Definition 7.3. We define the following sequence of co-ordinate transformations 

1 



(7.110) 
(7.111) 



3+ = r 2 w{l)(l - g- 1 / 2 a<r 1 )(l - q- l / 2 at-H) lf 



I -I 



3- 



iVr 1 )^ - g- 1/2 crf/)(i - q- l ' 2 at- l r l ) J—L : 

if i 



where the new variable f(t) is introduced. Clearly it is apparent that under the mapping I i— > l~ 1 we have 
the interchange 3+ o 3_ and that the product satisfies the relation (|7.27l) . i.e. 3+3- = 1X7=1 (1 — q^ 1 ^ 2 ajl ± ). 
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Let (T4 = q 2 s\. In addition we define the further variable p(t) so that 

r 1+ d n _ 1 2as 4 p- q~ n t~ q n a 2 slt- 1 
( ' P+ : ~gi/2[ n +|] t-t- 1 

ri_g» _ 1 2as 4 p - q n a 2 s\t - q- n t~ x 

Together these relations ensure that (I7.86|) is automatically satisfied. We further define <?(i) by 2p := g + g^ 1 ■ 
We shall find that /, /, g will be our primary variables. 

A consequence of these definitions and the previous Lemma is the following transformation formula. 

Corollary 7.5. Let us assume f ^= 0,oo, / ^ l^ 1 andw(f) ^= 0. There exists an invertible mapping between 
p(t) and X(t) with the form 

(/) 



(7-114) - 2s ^ +1 + ^ 2 =/^l-2/A- 

Proof. Employing the definitions (|7. 1 13[) . (I7.110|) and (|7.111[) in (|7.108[) we compute the relation 



(7.H5) 2s 4P = q -^a 3 -4[l + l {p _ 1){l _ f) ■ {p _ mf _ iy 

where we have shifted the arguments t H> g _1 ' 2 f of all variables. The relation (I7.114j) then follows by 
performing a partial fraction expansion of the above equation with respect to the variable l(t) (the former 
equation can be thought of as an expansion with respect to f(t)). We also have an alternative form of the 
above expression 

(7.H6) A = qV**L ± {g + g -i) _ 1 Mstg) f ^V)^ 

1 ' 2 CT4 2 Si {y g ' 2( 5 2 -l)/- s - 1 5 2( 5 2 -l)/- s -V 1 

□ 

Remark 7.6. What is surprising about this result in comparison to the situation with the g-linear lattice 
is that instead of a linear transformation between p and / we find one between p and A with / being an 
intermediary. Given 1 1 5[) it is clear there is no linear inversion possible for /. 

Our next task is to compute expressions for the four independent evaluations of the spectral matrix 
elements given on the left-hand sides of (|7. 104117. 1071) in terms of the variables introduced above. We will 
carry this out in stages of successive refinement. 

Lemma 7.2. The ratios of the evaluated spectral matrix elements satisfy the relations 
(7.117) 

a n 2as 4 p - g- n t - q n <x 2 s\t- x = W+(z 6 ; q~ 1/2 t) l+jz^g^H) = W-{El~ z^q-^H) z 6 ;g^ 2 t) 
a n 2as 4 p - q"a 2 s 2 t - q-'H' 1 %+{z 6 ;q- l / 2 t) 2B+(z 5 ; q l ' 2 t) 1+{El~ z 5 ;q- l / 2 t) W-{E 2+ z 6 ; qVH) ' 

Proof. Subtracting ([7T90]) from ([7792]) gives 

, 7ml /„+ jp- \ ( W-(E 2 u +z e ;qV 2 t) W + (z 5 ;gVH) \ 

whilst subtracting (|7.94[) from (|T.96[) gives 

y% + (z e ;q-i/ 2 t) ' ^ + (E 2 u -z 5 -q-^ 2 t) 



, 711Q , ,„+ jp- \ ( W + (z 6 ;q~i/ 2 t) W_(E 2 u -z 5 -q-VH) \ 
( 7 - 119 ) ri + \ E u x e ~ E u xs) = = ; . _ 1/2 ^ + ~ /r , 2 - rTTTTT P+- 



Taking ratios of these two gives 
r719n , r 1+ 1 + (z,-q^ 2 t)T + (E 2 u +z 6 ;q^ 2 t) 

(t.LZU) _ — ^ , _i/o.\t7- /r^2- 



ri_ % + (z 6] q-V 2 t)<Z + {E 2 u -z b -q-V 2 t) 

T + (z 6 ; g-^W^El-z,; q^'H) + l+jE^ z 5 ; q-^W+jz^ q^H) 
^ + {z^qy 2 t)W-{El + Z Q] q l / 2 t)+^ + {E 2 u + z Q -q i m)W + {z 5 -q^H) 
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Likewise subtracting (|7.91[) from 93[) yields 

rri-m (w+ r- \ f ^ + (E 2 u + ^,g 1/2 t) , 1 + (z 5 ;g^t) \ 

and (fTM) from ([71371 

( 7 - 122 ) r i- ( E u x e -E u x 5 ) = -[ — — — -pr- + — ,_ — — — p-. 



Taking the ratio of these latter two equations gives 

(7 1231 ri+ - ^+( z e;q- 1/2 t)W-(El-z 5 ;g- 1 /H) 
' ' n- W + (z 5 ;qVH)W-(El+z 6 ;qy 2 t) 

X+(z 5 ; q^W-iEl+Ze; q^H) + T+(E 2 u +z 6 ; qVH)W+(z B ; q^H) 



1 + {z 6 ; q-VH)V3^(El-z 5 ; q-VH) + 1+(E'tz 5 : q-V 2 t)W+{z 6 ; q~V 2 t) ' 



The factors appearing on the right-hand sides of these two ratios constructed have simple evaluations, the 

first factor being 

(7.124) 

% + {z 5 -q 1 / 2 t)1+{E 2 +z^q 1 l 2 t) _ a 2 n (q 2 t 2 - a 2 ){a 2 t 2 - l){at - l){at - t^jgt - at^qt - al) 
T + (^ 6 ;g- 1/2 t)1 + (£;2"%;g- 1/2 t) ~ < (t 2 - a 2 )(a 2 t 2 - q 2 )(at - ql){at - qh l ){t - ah l ){t - al) ' 

Employing this evaluation in the first ratio (|7.120[) yields the evaluation for the second factor 
T + (z 6 ; q-V 2 t)QB-(El-z 5 ; q~ x l 2 t) + 1 + (El~ zs;q-VH)W + {z*\q-VH) 



(7.125) 



1+(z 5 ; q^ 2 t)W-(E 2 +z 6 ; q V 2 t) + 1+(E 2 +z 6 ; qV 2 t)W + (z 5 ; qVH) 

_ a n (t 2 - a 2 )(a 2 t 2 - q 2 ) {at - ql){at - ql^){t - at 1 )^ - al) 2as 4 p - q~ n t - q n a 2 s 2 4 t- 1 



i n (q 2 t 2 - a 2 ){a 2 t 2 - 1) ( a f - l)(at - h^qt - ah^qt - at) 2as 4 p - q n a 2 s 2 t - q^H- 1 ' 



Employing this latter factor in the second ratio (|7.123[) we can construct two relations involving perfect 
squares, upon using (|7.103[) . After taking the square roots we need to resolve the sign ambiguity. The final 
result is then (|7.117|) . □ 

Remark 7.7. In fact we can make separate evaluations of the numerator and denominator of the last ratio, 
which gives us 

(7.126) 1+(z 5 ;q^ 2 t)W-(El + z e ; q^t) + 1+(E 2 u +z 6 ; q^ 2 t)W+{z 5 ; q^ 2 t) 
2 a \n + -1 

= " L . 21 (q 2 t 2 -a 2 )(a 2 t 2 -l)(a 2 -q)(at-i)(at-r 1 )(qt~al)(qt-ar 1 ) \2as 4 p - q n a 2 s\t - g - "* -1 ! , 
and 

(7.127) T + (z 6 ; q-VHytB-ffi- z 5 ; q~ x l 2 t) + X+(^- Z5 ; q -VH)*B+(z«\ q^'H) 
2d In + -] 

= " 1 _ . 2 \ t 2 -a 2 )(aH 2 -q 2 )(a 2 -q) (at - ql ) (at - ql~ x ) (t - al) (t - al' 1 ) \2as 4 p - q~ n t - q n a 2 s 2 A- 1 ] . 
q i a°t i 

In order to proceed any further we will require representations of the individual evaluated spectral coef- 
ficients which we give as two distinct partial-fraction expansions. We offer these without proofs as they are 
the outcome of straightforward although lengthy computations. 
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Lemma 7.3. The evaluated spectral coefficients possess the rational function forms in the /, I variables 



{r . m 2([n][n+ . R 2Mg^| 

= -S-'rV - afKl'atn - <,-/) - + °T l mf-\<*t - «-"/)(«"<>* - /)(! + H) 

(.J v(i ' / 

+ a-\r 2 (q n at - /) [ - <T"a 6 (/ 2 - l)(g"ai + /) - q-^ataeaj 

+ cr n - 5 / 2 a 2 a 3 f 2 + (q- n+2 + q n - 2 a 4 a 4 )f(-q- n - 2 + q- 4 ataj)] , 



(7,29) «[„][„+ i^^gi^J 

= _ -i >-2 ft - <*/)(<** - qf){q n ^ - q~ n tf) w(f) 
(a-tf) ~ ' 
+ a- 1 6r 6 /- 1 (at - g- n+1 /)(« n " 1 a* - /)(* + 
_ [n](g-t 2 )i 3 (f-a/) m(ar') 

9 (a-*/) (a-tZ)(a- tf- 1 ) 
+ a" 1 *" 1 / - V _1 a* - /) [ - <T n tcr 6 (/ 2 - l)( g "trf + «/) - g" n (g - t 2 )a 6 f(q n t + af) 

+ q- 2 - 2n t 2 f(-q 2 t + «" W) - g-^Ww/ + ^""a 2 *^/ 2 - - ^W)] , 



(,130) «[„][„+ i^f^i^ 

= q-Wf-^at f)(q n+1 / 2 ta 6 9~ n ~ 1/2 af) J"^ - a^aj-^at q- n f)(q n at f){l + H) 

+ /" V«* - /) [q- 3 - n (q n - l)a<7 4 /(l + at/) + q- 1/2 ta 6 a!f - q- n ~ 5 / 2 aa s f 2 

{q n+1 ' 2 ta, q-o-Wafrq-"- 1 ' 2 + q^aj 2 )} , 



(7 131) 2t\n]\n I Mfi 

(7.131) 2 *N[n+ 2 ja n ^ +( ^_^^_ 1/2 ^ 

= _ 1/2 f -_2 ft - «/)(crt - g/)(g" +1 / 2 a 6 - q- n -V 2 atf) w(f) 
C1 (q-atf) (/_f)(/_f-i) 

- a-'aj-^at - q~ n+1 f^q^at - f)(l + t 1 ) 
[n](q - t 2 )aH 3 (at - qf) w(qa~ 1 t- 1 ) 

q 2 (q - atf) (q - atl)(q - ath 1 ) 

+ i" 1 / - V" 1 "* - f) [q^a-Haeif 2 - l^q^at + f) + q-*~ n (q - t 2 )<jj(q n a 3 t + q 2 f) 

+ q-y 2 t 2 a 6 aj - q- n ~ 3 / 2 ata 3 f 2 - q- 4 - 2n a(qt 2 + q 2n a A )f{-q 2 t + q n aa i f)] . 

In addition the evaluations at the advanced co-ordinate have alternative forms. 

48 



Lemma 7.4. The advanced evaluated spectral coefficients are rational functions in the g, f variables 
W+iz^q^t) 



(7.132) 2t[n][n+ i]a n 



_n-i (g" tos 4.9 1 ^ !) w(s 4 1 g) _„_! (q n tas 4 g - if w(s 4 l g v ) 



g-"- 1 a(.9 2 + .g- 



.9 -.9 



5-5" 



/-S4 1 <T 



q- 2 "- 6 ^- 1 [q 6 + t 2 q 4n a 6 a 2 - t 2 q 3n+1 a 4 a 4 (q 2 + o 4 



qV 2 aa 3 + ta 4 (q 4 2q n + 3 a 2 + q 2n a*a 4 )\ (g + g" 1 ) 
2 +a 4 )+g 2 "+ 2 a 2 (tV+a> 4 
<f + 4 a 2 (2< Z + a 2 ) + q 2n+5/2 a 3 t(-qa 3 + a x a 4 ) 



(7.133) 2i[ra][ra- 



gu_(^ 6 ;g 1/2 

2J a ™ rr ,-7^2+ 



5 ™ 1 as 4 



(g^ias^ -1 - l) 2 w(s4 _1 5) 



+ q n 1 a(g 2 +g ~)-q 



2^ _„-„-15/2 aS4 -l 



-n-i (q n tas 4 g - 1) w(s 4 *</ *) 
5-.9- 1 f-s 4 ^g-i 
taa 4 (q 9 > 2 - 2q n+9 / 2 + g 2 "+ 5 / 2 a 4 )j (g + g' 1 ) 



. q -2n-U/2 a 



-q 9 ' 2 - q in+1 l 2 a 2 t 2 a 2 + q 3n+1 / 2 a 2 t 2 a 4 (q 2 + a 4 ) - q 2n + 5 / 2 (l + a 2 t 2 )o 4 



+q n+7 / 2 (2q + °2) - q 2n+2 at(-qa 3 + a,a 4 ) 



At this stage we have accumulated enough results to deduce the dynamical equations for our system on 
the deformation lattice. 

Proposition 7.9. The t-evolution of our system in the variables p, A, / is given by a pair of coupled first 
order equations the first of which is 



(7.134) 



2qs i p-t 2 f- 1 -q 2 t- 2 sif = q(t-af)(qf -at) 
2s 4 p - f- 1 - s 2 J (tf -a)(q- atf) 

aH 2 (q - t 2 ) (qf - at) /^(ga" 1 *- 1 ) 



+ 



(q - a 2 ) (q - atf) qatw(f) - (qf - at)(q - atf)f[2s 4 p - /-I - s 2 J] 



qt 2 (q-t 2 )(t-af) 



f 2 w(at-i) 



(q - a 2 ) (tf - a) atw(f) - (tf - a)(t - af)f[2s 4 p - /-l - s 2 f] ' 
The auxiliary equation for the leading coefficient of the polynomial is 



(7.135) 



t 



2ln 

7 2 
in 



2 t 2 + q 2n a 2 s\ - 2q n ats 4 p 



2 l + q 2n a 2 t 2 s 2 4 -2q n ats 4 p 



Proof. The first of these equations is derived utilising the following steps. We substitute the expression for 
ri+a n /p+ given by (|7.112p and the transformation formulae for |± as given in (|7.110|7.111|) into the solution 
at the retarded co-ordinate (|7.109j) . This yields an equation involving p on one hand and /, / on the other 
hand. Now we perform a partial fraction expansion of this with respect to I and this produces the following 
expression 



(7.136) 2s 4 p = q-H 2 f- 1 +qt~ 2 s 2 4 f 



(t-af)(qf -at) w(f) 



f(tf-a)(q-atf) p + l-2f\ 
aH 2 (q-t 2 )(qf - at) w(qa- 1 t- 1 ) t 2 (q - t 2 )(t - af) w(at- 1 ) 



q(q- a 2 )(q- atf) q 2 + a 2 t 2 - 2qat\ (q - a 2 )(tf - a) a 2 + t 2 - 2atX 
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which is a simple function of A. Then one substitutes for A using the inversion of the transformation (|7.114[) 
at the retarded time and the result is (|7.134l) . Alternatively one can prove this formula by substituting 
CLUl and dUnH) into (|TTT27|) . 

We employ the variable transformation (|7.112l) and note that the left-hand side of this expression has 
been evaluated in (|7.81[> and (|7.83l) . Equating these two forms gives (|7.135p . □ 

Remark 7.8. The right-hand side of formula (|7.134[) exhibits apparent poles at a = tf and q = atf. However 
this is not the case as the former is cancelled by opposing contributions from the first and third terms, while 
the latter is cancelled by contributions from the first and second term. 

Clearly (|7.134p is not manifestly invertible for p, however it is possible to construct a linear equation for 
the retarded variables by switching from p to A. 

Proposition 7.10. The inverse to \1.131$ is given as a relation for X, which we give in two alternative 
forms 



(7.137) A 



(/ - q n at)(t - q n as 2 J) 



(tf - q n a)(l - q n ats 2 J) 



t 4 + f 2 



2t 2 f 



t 2 (t 2 - \)s\ 
~2q n + 1 / 2 [n + i] 



(f-q n at) q An a i w(q- n a- 1 t- 1 s± 2 ) - Z (t - q n as 2 4 f) 



w(q n at x ) 



' 4 (1 - q n ats 2 J) 1 + q 2n a 2 t 2 s 2 - 2q n ats 4 p 4 (tf _ q n a ) t 2 + q 2n a 2 s\ - 2q n ats 4 p 



or 



(7.138) 



(g - q n ats 4 )(t - q n as 4 g) w(s 4 1 g) g 2 (l - q n ats 4 g)(tg - q n as 4 ) w(s 4 1 g x ) 



2(g 2 - 1)(1 - q n ats 4 g)(tg - q n as 4 ) f - s j l g 2(g 2 - l)(g - q n ats 4 )(t - q n as 4 g) f - 



2at 3 



a 2 t(t 2 - 1)4 
2q 1 / 2 [n + i] 



(74 



q 2n as 4 



g V2 ato3 + g5/2 (t 2_ 1){n+ l } 

w(q- n g- 1 t- 1 s 4 2 ) 
(1 — q n ats 4 g)(l — q n ats 4 g~ 1 ) 



2t 2 



.s 4 



\g + g~ 



w(q n at 1 ) 



(t — q n as 4 g)(t — q n as 4 g 1 ) 



Proof. One can solve the compatibility relation (|6.67[) for the (1,2) component of A* n (z;q x l 2 t) and find 
% + {z;q- l IH) 



(7.139) 1 + (z;q-^ 2 t) 



X(z, t) W 2 (R + AvS) (q l l 2 z; t)(R - AvS) (q l l 2 z; t) 

9M<z 1/2 ^-(^ 1/2 ^)^; 9 17 ^ ^ 

+D\-(q 1/2 z; t)^ + (q-^ 2 z; t)W+(z; q 1/2 t) - ^ + (q^ 2 z; t)«R_ (q- 1/2 z; t)2H_ (z; q 1/2 t) 

To simplify the calculations of the spectral matrix elements we employ an alternative parameterisation to 
that of ([737]) . ([73511 . ([7351 



(7.140) W+(z;q 1/2 t) 



at(z 2 ~ l)(qt ~ az)(q 2n+ Ha 6 z - a) 
' z 2 (l-a 2 t 2 )(q-a 2 )(q 2n + 1 t 2 d 6 - 1) 
q 2 at(z 2 - l)(z - at)(q 2n ata 6 z - 1) 
z 2 (q 2 t 2 - a 2 )(q - a 2 )(q 2n +H 2 a (i - 1) 

(z + l)(z - cd)(qt - az)(q 2 - q 2n a 2 t 2 a 4 z) 
+ 2qtz 2 (q 2 ~ q 2n a 2 t 2 a 4 ) 

(z - l)(z - at)(qt - az)(q 2 + q 2n a 2 t 2 a 4 z) 



2B+(crf;q 1/2 i) 
W^q- 1 ^- 1 ^ 1 ^) 
w(l) 



2qtz 2 (q 2 - q 2n a 2 t 2 a 4 ) 



w(-l) 



-n-3f Jl 



(z — l)(z — at)(qt — az) 
at 2 z 3 (q 2 — q 2n ct 2 t 2 o- 4 ) 



+ q 4n a 4 t 4 a 2 z 2 + q 2n+3 / 2 a 2 t 2 (qa 3 z + a x a 4 z - q x / 2 (l + z 2 )) 



5(1 



, _ . , ^ / . 1/2.N Qit(z 2 - l)(qtz - a)(q 2n+1 ta 6 - az) 1/2 
(7.141) g t) = ^ _ aH2)(q _ a2){q2n+H2 - 6 _ 1} W + (at; q t) 

q 2 at(z 2 — 1)(1 - atz)(z — q 2n ata e ) 
+ z 2 (q 2 t 2 - a 2 ){q - a 2 )(q 2n +H 2 d 6 - 1) 

(z + 1)(1 - atz)(qtz - a)(q 2 z - q 2n a 2 t 2 o±) 



2qtz 2 (q 2 - q 2n a 2 t 2 a 4 ) 
(z - 1)(1 - atz)(qtz - a)(q 2 z + q 2n a 2 t 2 a 4 ) 



W-iq^ar 1 ^ 1 ^) 
w(l) 



2qtz 2 (q 2 - q 2n a 2 t 2 a A ) 



w(-l) 



-n-3( J2 



and 

(7.142) l+frq^H) 



(z 2 — 1)(1 — atz){qtz — a) 
at 2 z 3 (q 2 — q 2n a 2 t 2 GA) 

at(z - z- 1 ) 



1„2 , „4n„4 + 4„2 , „2n+3/2„ 2 + 2 ( „„ „ , 

z + q a t <7 4 + q 'at [qcr^z + 



o~io~4,z — q 



l ' 2 {l + z 2 )) 



(q~a 2 )(qt 2 -l) 
(qtz — a)(qtz _1 



P+ 



W+(at;q^ 2 t) 



(a 2 t 2 - 1) n,_|(«- 1 /2 C rf + g i/a a -x t -i) + rQ _ 

2 (z - at)^- 1 - at) W-iq^ab^iq^H) 



+q 



(a 2 - q 2 t 2 ) n,_i(g-V2 a t-i + qWcrH) + r 0l - 



(7.143) 1_(z; 9 1 / 2 <) = 



at(z — z - 1 ) 1 
{q-a 2 )(qt 2 - 1) p+ 



( ^~ ( ^ij~ a) (a-^(g- 1/2 °* + ^a" 1 *" 1 ) + oar-Cat; ? 1/2 t) 

(cr — qH') 

This particular parameterisation implies that the numerator of the right-hand side of (|7.139p manifestly 
contains a factor of (z — z~ 1 )(qtz — a){qtz~ l — a)(z — ai)(z _1 — at) which is present in the denominator. 
This ensures that the ratio is linear in x, as it must. Finding the zero of this linear polynomial then gives A, 
which after further substantial manipulation and simplification yields (|7. 138|) and (I7.137[) . □ 

A crucial fact enabling further progress is the following factorisation formula for a quantity that will 
subsequently figure prominently in certain discriminants. 

Lemma 7.5. The bi-quadratic in p, A 

(7.144) 16cr 4 (p 2 A 2 - p 2 - A 2 ) - 8s 4 {q 2 + qa 2 + cr 4 )pA + 8q 1/2 s i (qa 1 + a 3 )p 

+ 8q- 1/2 (<j 1 a 4 + <?cr 3 )A + (q- o 2 ) 2 - 4cricr 3 + 2cr 4 - 2g~ 1 cr 2 cr 4 + <T 2 er 2 , 
is a perfect square which can be given in either of two ways. In the first way this is the square of 
(7.145) 

q?f-*-q L ' i 1?a l r l + q- l 'H- i v 3 f-q-H-*v 4 p + ^ 

(1 - s 2 J 2 )(f - q n at)(t - q n as 2 f) w(f) 



2qs i (t- 2 f-t 2 f- 1 )p-q 
q 1 ' 2 at{t 2 - l)s\ 



f 2 (tf - q"a)(l - q n ats 2 J) 1 + s\f 2 - 2s 4 fp 
g 2 "a 2 sl(l - q 2n a 2 t 2 sl)(f - q n at)w(q- n a- 1 t- 1 s^ 2 ) 
(1 - q n ats 2 J)(l + q 2n a 2 t 2 s\ - 2q n ats±p) 

q- 2n a- 2 s- 3 (t 2 - q 2n a 2 sl)(t - q n as 2 J)w(q n at- 1 ) 
(tf - q n a)(t 2 + q 2n a 2 s\ - 2q n ats 4 p) 
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and in the second way this is the square of 



(7.146) gVVat- 2 / - g^Vx* 3 /" 1 - {t 2 1 2 + gV 2 * 2 ) [f - 1 + q~\q + a 2 )^ 1 *"^ - t 2 )f] 

q(f 2 - l)(qf - at)(t - af) w(f) 



+ 2A(t 2 /- 1 - gV 2 * 2 ,/) + 
t(g-< 2 ) 



qa{q — a 2 ) 



f 2 {tf-a)(q-atf) p + l-2f\ 
t - af q 2 (t 2 - a 2 )w(at- 1 ) qf - at a 4 (q 2 - a 2 t 2 )w(qa- 1 t- 1 ) 



tf -a a 2 +t 2 - 2atX q - atf q 2 + a 2 t 2 - 2qatX 



Proof. In the first way we substitute (17. 138[) into (|7.144l) . whereas in the second way we substitute (|7. 136j) 
into (|7.144[) . After considerable simplification we arrive at the two results. □ 



We are now in a position to derive the evolution equations for the /-variable, firstly in the advanced 
direction. 

Proposition 7.11. Assuming f ^ 0, tf — a ^ 0, and q — atf ^ we have the forward evolution for f 



(7.147) / = 



2t 2 {t 2 - l)s 2 

n-l/2I n+ U 



q 



q An a i s\w{q- n a- 1 t~ 1 sl 2 ) 



w(q n at 1 ) 



1 + q 2n a 2 t 2 s\ - 2q n ats 4 p t 2 + q 2n a 2 s\ ~ 2q n ats 4 p 
-at- q-\l - 2t 2 )<7 4 - q(l - 2t- 2 ) + 2A [2qs 4P - i 2 /" 1 + q 2 r 2 s 2 f] 
q 1/2 a 3 t- 2 f + q^a^r 1 + (t 2 f~ 2 + f^^l) [P -l + <f\q + a 2 )^ 1 ^ - t 2 )f] 
q(f 2 ~ l)(qf - at)(t - af) w(f) 



t(q-t 2 ) 
qa{q — a 2 ) 



f 2 (tf-a)(q-atf) f 2 + 1 - 2/A 

t - af q 2 (t 2 - a 2 )w(at _1 ) qf - at a 4 (q 2 - a 2 t 2 )w{qa~ 1 r l ) 
tf -a a 2 +t 2 - 2atX q - atf q 2 + a 2 t 2 - 2qatX 



' 2q 1 ' 2 t{t 2 - I) 



q 2n a i slw{q- n a- 1 t- 1 sp) 



as\w(q n at r ) 



1 + q 2n a 2 t 2 s\ - 2q n ats 4 p t 2 + q 2n a 2 s\ - 2q n ats 4 p 
2q^ 2 a~H- 3 {n + \}{t 2 - 1) - 2q-H- 2 (q 1 / 2 a 3 - 2t 2 a 4 X) + Aqr 2 s 4 p 



where p is given by {7.13$ - 



Proof. Whilst (|7.138l) can be primarily viewed as a linear equation for A it can also be viewed as a quadratic 
equation for /, and as such, possesses a discriminant which contains, apart from explicit squared factors, the 
factor (I7.144[) . By substituting for p using (|7.136l) we can employ the result (I7.146[) . and effect a factorisation 
of the quadratic into linear factors. The choice of the factors can be settled by consideration of the known 
solution for n — and it transpires that the negative branch is appropriate. This yields l|7.147l) . □ 



Our last task is to derive the inverse to (|7.147[) and complete the system of evolution equations. 
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Proposition 7.12. The inverse of the evolution for f is given by 



(7.148) / 



2t 2 (q-t 2 ) 



wiat- 1 ) 



a 2 +t 2 



2at\ 



q-'a 4 - 



w(qa t 1 ) 
+ a 2 t 2 - 2qatX 



ASipX + 2s A {t- 2 f - t 2 f- X )p - t 2 f- 2 + q- 1 ' 2 t 2 o- 1 f- 1 - q- X G 2 - q 
-2 4 



3/2 



*~ W + q- 2 t- 2 a 4 f 2 



q 2 r 2 (-2q + t 2 - t 4 )<7 4 - q^t^^l + t 2 ) - 2t 4 ) - gV a {n + ^a" 1 *" 3 ^ 2 - 1)(/ + t 4 /" *) 



(1 - s 2 J 2 )(f - q n at)(t - q n as 2 f) 



2i 

w(f) 



at(t 2 



f 2 (tf-q"a)(l-q n ats 2 f) 
q 2n a 2 s\{\ - q 



1 + sif 2 - 2s 4 fp 



-In, 



2n a 2 s\{\ - q 2n a 2 t 2 s 2 4 )(f - q n at)w{q- n a-H- 1 s 4 2 ) 



(1 - q n atsjf)(l + q 2n aH 2 s\ - 2q n ats 4 p) 



- 3 (t 2 - q 2n a 2 s 2 )(t - q n as 2 4 f)w(q n at~ 1 ) 



(tf - q n a){t 2 + q 2n a 2 s\ 



j 2at(q - t 2 ) 
| q — a 2 



w(at 1 ) 



2q n ats 4 p) 
a 2 w(qa~ 1 t~ 1 ) 



- 1 2 - 2at\ q 2 + a 2 t 2 - 2qatX 
2q- 2 a- 1 r :i (q - t 2 )(q + a 2 )a 4 + 4s 4 p + 4g" W 4 A 



where X is given by \1.138\ . 

Proof. In a similar manner to the previous proof we can solve (|7.136j) for /. On appearance this polynomial 
should be a quintic in /, however in line with Remark l7.8l we observe this contains a factor of f(tf—a)(q—atf) 
and we only have to deal with a quadratic. Upon examining this quadratic we find that the discriminant 
contains, again, the factor (|7.144l) as the only manifestly square-free factor. Upon substituting for A using 
(|7.138l) we can now employ (|7.145|) and factorise the quadratic. As in the previous case consideration of the 
n = solution resolves the sign ambiguity in favour of the negative root, yielding (|7. 148|) . □ 

7.4. Seed Solution. We can now make contact with the earlier theory characterising the moments, as given 
in Subsection 17.11 through a study of the classical "seed" solutions to the coupled recurrence system given 
above. We now append a subscript on the variables to indicate the n— value. 

Proposition 7.13. The recurrence relations \1.131$ and admit the classical "seed" solution at n = 



(7.149) 
and 

(7.150) 



hit) 



-1/2 



at, 



2qs 4 p (t) 



C(q 1 / 2 a- 1 t + o-iq-^at-^mit) + D(g 3 / 2 a _1 i _1 + a iq - s / 2 at)m{q-H) 



Cm(t) + Dm(q-H) 
where m(t) is the general solution to the second-order q-difference equation 

(7.151) (t 2 - q^wiat'^miq 1 / 2 ^ 

+ 5 [{t 2 - q^wiat' 1 ) + q 3 (t 2 - ljwiq^at) - q- 3 (q - a 2 ){q 2 - a 4 a 2 )(t 2 - l)(t 2 - q){t 2 - q 2 )} m{q- 1/2 t) 

c 

+ 7?t^ 3 (* 2 - iM<rW)m(<r 3/2 i) = o. 

Proof. It is clear that both (|7.130[) and (|7.131l) both vanish when (|7. 149|) holds at n = 0, therefore together 
they satisfy (|7.117[) . Given (|7.149l) we observe that two terms vanish in (|7. 134|) and po(t) satisfies the discrete 
Riccati equation 



(7.152) 



2s 4 po — a t 



.-1 2 

at s A 



qt^wiat- 1 ) [2qs 4 p - q 2 a - 1 t- 1 



ats 2 ~\ 



at(q — t 2 )(q — a 2 ) [2qs 4 po — q 2 a 1 t 1 — ats 2 ] + q 4 w(q 1 at) 
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Alternatively we can derive this relation by specialising (|7.137[) under (|7.149p which yields 

7.153 t +a -2at\ = — — ^ — — — , 

and then using (|7.114l) to substitute for Ao- Making the standard linearising transformation 

Am(t) + Bm(q~H) 

we find, in our solution for the decoupling factors, that A = (qa~ 1 t + d4q~ 1 at~ 1 )C and B = [q 2 arH~ l + 
o^q~ 2 a£)D and arrive at (|7.15ip . □ 

Remark 7.9. Our explicit initial orthogonal polynomial variables are one specialisation of the classical so- 
lutions given above. As remarked in the proof of Proposition 17.31 we have fio = V(\q) and vq = W{\q) at 
n — 0, where Ao is the unique zero of the right-hand side of (|7.18l) 

(7.155) 2[±]A = -g-V» ffl + q^o, - [l] m °' + + m °'- . 

m ,o 

We easily see that 3±,o(*) = ^ 3 IIj=i(l — 9~ 1/ ' 2a j'o tl ) an d from the transformations (17. 11017. lllj) we deduce 
that /o(t) = q^ x ^ 2 at. The elements of the spectral matrix Aq are 

(7.156) %8 +fi (z; f) = z~ 3 w(z)(l - q- 1/2 atz){\ - q'^at^z) = W + AyV, 

AyV, 



(7.157) 


2B_. (z;f) = z 3 u)(z~ 1 )(l - g- 1 / 2 a te" 1 )(l - q 




(7.158) 


T +l0 (z;t) = 2a [\]{z- Z - l ){z- 


-/o)(l-io X 2 




(7.159) 


T_,o(«;t)=0, 






whilst those of the deformation matrix are 






(7.160) 


9t+,o(^;t) = _ . (1 
ao7o 


<T 1/2 atz)(l - 


-q-^atz- 1 ), 


(7.161) 


9t-,o(^;*) = - ~ (1 

ao7o 


q- 1 / 2 at- 1 z)(l- (? - 1 / 2 a t- 1 z- 1 ), 


(7.162) 


qj+.ofo*) =2 (? - 1 / 2 aa 


. 7o 7o_ 


5 


(7.163) 


«P_,o(«;t) = 0. 







We also note from ([7335]) and (|3.16|3.10|l that 

_ ,.v (1 + <7" 1 a 2 s^ 2 )rn ,o(g _1 - (1 + 9a 2 si^ 2 )mo,o(i) 
(7.164) 2as 4 p {t) = — — - yj^— — , 

which is precisely the case of (|7. 15Q[) with D = —q~ 1 t 2 C ' . With this solution we also deduce that the moment 
recurrence (I7.78[) for m 0j o coincides with (|7.151l) . 



At the beginning of this section, in Remark 17.21 we noted that some explicit solutions to the moment 
recurrences were given in [97j . This knowledge enables the possibility of checking any aspect of the foregoing 
theory either symbolically or numerically to essentially unlimited precision. We wish to report that extensive 
checking of all the key relations for the cases n — and n = 1 has been carried out employing Mathematica 
code and utilities where possible by exact symbolic means or if not by numerical means to a level of one part 
in 10 20 — > 10 30 or better at random exact values of the input parameters. These checks have also utilised a 
body of unreported work [96] covering the system of Laguerre-Freud recurrences under n ^ n + 1 for general 
n e Z>o- 

We now come to the question regarding the identity of this system in the Sakai scheme. Whilst not 
providing an explicit transformation from our parameters and variables to those appearing in the canonical 
coupled ^-difference system [33l [86] , we offer unambiguous evidence that it is one of the classical yet full 
parameter set cases of the £7 g-Painleve system as can be seen from the following inspection of their 
classical solutions [7IJ [51] [52], and the recent systematic study of "two-Casorati" determinantal forms of 
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their classical solutions [69 . Clearly the r- functions of our construction, given by (|3.10p and p. lip , are of 
the "two-Casorati" determinant form because if we choose b — a r , a — a s with r, s € {1, 2, 3, 4} say and this 
is employed in the definition (|3.9[) we have a moment determinant with elements to Ui o(<7 j Q>r, q k o,s), j, k £ Z>o 
where mo^ given by any (/-constant linear combination of the two gW-j solutions appearing in (I7.9[) . 

Several tasks, that have arisen in the course of our study, remain unfinished and we conclude by de- 
tailing them and the issues involved. Clarification of the explicit relationship of the evolution system 
(|7. 13417. 13717. 14717. 148[) with the canonical g-difference equations as say given in (5TJ [52] is required. It 
remains to complete the construction of the n — > n + 1 recurrence relations which occupy a special place 
in our approach, but figure as Schlesinger transformations in the integrable theory. In fact our framework 
can easily treat the Schlesinger transformations of this example (or any for that matter), and in particular 
the n — > n + 1 recurrences, or the aj, a i— > qa,j, qa, j = 1, . . . , 4 mappings. These latter transformations are 
manifested in our context as specialised Christoffel-Uvarov transformations, however we will postpone this 
undertaking. 

Significant progress has been made in finding the analog of an isomonodromic system for the elliptic 
Painleve courtesy of a preprint by Eric Rains |85j , and in the work by Yamada |98j . The approach taken in 
this former work is very much in the spirit of the present study and it would be natural to expect that a 
limiting case of the results reported there would correspond our own. Very recently Yamada [99] has given 
Lax pairs for the E^\ Ey and Eg q-Painleve equations by reformulating the E^ elliptic Painleve Lax 
pair system and taking limits -Eg 1 '' — ¥ Ej 1 ^ — > E^ . One would expect expect that a gauge transformation 
and co-ordinate transformations would link our Lax Pair with his. 
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